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ABSTRACT 

We examine the physical properties and evolutionary stages of a sample of 17 clumps within 8 
Infrared Dark Clouds (IRDCs) by combining existing infrared, millimeter, and radio data with new 
Bolocam Galactic Plane Survey (BGPS) 1.1 mm data, VLA radio continuum data, and HHT dense 
gas (HCO + and N2H" 1 ") spectroscopic data. We combine literature studies of star formation tracers 
and dust temperatures within IRDCs with our search for ultra-compact (UC) H II regions to discuss a 
possible evolutionary sequence for IRDC clumps. In addition, we perform an analysis of mass tracers 
in IRDCs and find that 8 /im extinction masses and 1.1 mm Bolocam Galactic Plane Survey (BGPS) 
masses are complementary mass tracers in IRDCs except for the most active clumps (notably those 
containing UCH II regions), for which both mass tracers suffer biases. We find that the measured 
virial masses in IRDC clumps are uniformly higher than the measured dust continuum masses on the 
scale of ~ 1 pc. We use 13 CO, HCO + , and N2H" 1 " to study the molecular gas properties of IRDCs and 
do not see any evidence of chemical differentiation between hot and cold clumps on the scale of ~ 1 
pc. However, both HCO + and N2H" 1 " are brighter in active clumps, due to an increase in temperature 
and/or density. We report the identification of four UCH II regions embedded within IRDC clumps 
and find that UCH II regions are associated with bright (> 1 Jy) 24 /im point sources, and that 
the brightest UCH II regions are associated with "diffuse red clumps" (an extended enhancement 
at 8 /im). The broad stages of the discussed evolutionary sequence (from a quiescent clump to an 
embedded H II region) are supported by literature dust temperature estimates; however, no sequential 
nature can be inferred between the individual star formation tracers. 
Subject headings: dust, extinction - ISM: clouds - stars: formation and pre-main sequence 



1. INTRODUCTION 

Despite their profound impact on the galactic envi- 
ronment and enrichment, the protostellar evolution (and 
birth) of massive stars (M > 8 Mq) remains elusive. Di- 
rect observations of the earliest evolutionary stages of 
massive stars is a challenging task. Massive stars evolve 
more quickly, are more rare, and consequently farther 
away on average than their low-mass counterparts. In 
addition, massive stars, once formed, quickly heat and 
ionize their environment, disrupting their natal molecu- 
lar cloud. Much like their low-mass counterparts, mas- 
sive stars are thought to form in cold, dense molecu- 
lar clouds. Massive stars seem to form predomin antly 
in clusters (|Lada fc Ladal 120031 : IdeWit et all 120051 ) and 
therefore, their natal molecular cloud should be more 
massive and denser than the natal clouds of isolated 
low-mass stars. Infrared Dark Clouds (IRDCs) have 
been proposed to b e the progenitors of massive stars and 
their host clusters (lEpn et al.l [l99l iCarev et all 119981 : 
iRathborne et al.ll2006jf) ~ 

IRDCs were identified as dark extinction features 
against the bright mid-infrared gal actic plane with 
the MSX and ISO space satellites ((Egan et all 119981: 
iPerault et ailll99l lOmont et al 1 [200l. IRDCs were cat- 
alogued from the MSX 8 /xm data bv lSimon et all ((2006a) 
and from the GLIMPS E 8 /tm data more reccently by 
IPeretto fc Fullerl ((2009). The strong continuum extinc- 
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tion and lack of infrared emission by these clouds belie 
their cold temperatures and high column densitites. Pre- 
vious studies have confirmed that these cloud s are indeed 
cold ( T < 25 K) and dense (n H > 10 5 cm" 3 ) (|Egan et all 
Il998t ICarev et al.l 119981 12000ft . More recently, studies 
have found kinematic distances to these clouds and de- 
tcrmi ned that their masses li e in the range of 10 2 to 10 4 
M Q (|Rathborne et alj 12006ft . iRagan et all (|2009ft find a 
clump mass spectra (over the range of clump masses from 
30 to 3000 M ) in IRDCs with a slope of about 1.76 ± 
0.05. 

There is compelling evidence to suggest that IRDCs 
are the cold, dense precursors to stellar clusters, and 
thus massive stars. Millimeter continuum maps show ex- 
tended cold dust emission, matching the morphology of 
the IRDCs and surrounding compact clumpJI- These 
clumps fall into two categories: active and quiescent. 
The active clumps are characterized by some or all of 
the following: 1) lu minous, embedded 24 /im emission, 
2) "green fuzzies," ((Chambers et alj 12009ft an extended 
4.5 /im enhancement believed to be caused by shocked 
H2, thus indicative of outflows (also known as Extended 
Green Objects, EGOs lCvganowski et aUl2008D , 3) maser 
emission, and 4) ultra-compact (UC) H II regions. A qui- 
escent clump is characterized as a cold, dense millimeter 
peak with none of the aforementioned signs of active star 

1 In this paper, we use the term "core" to describe a small, dense 
object within an IRDC clump, with a size of order 0.05 pc which 
will likely form a single star or stellar system, "clump" to describe 
a dense object within an IRDC with a size of order ~ 1 pc and 
a mass 10 2 - 10 3 Mq, and "cloud" for an IRDC which has a size 
scale of order 10 pc, and a mass of order 10 3 - 10 4 Mq. 
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formation. It is believed that the active clumps are cur- 
rently forming stars while the quiescent clumps are at an 
earlier stage of star formation: the cold, dense precursor 
before the star has "turned on." 

IRDCs are the densest clumps and fila ments embed- 
ded wi thin Giant Molecular Clouds (GMCs: ISimon et al.1 
2006b) with the favorable viewing condition of being on 
the near-side of a bright mid-IR background. They are 
preferentially found toward the Galaxy's largest star- 
forming region, the 5 kpc ring. IRDCs have typical sizes 
D ~ 1-10 pc, densities nff > 10 4 cm " 3 , and masses M 
- 10 2 -10 4 M (jRathborne et al.ll20"oll . IRDCs are dense 
condesations embedded within GMCs (D ~ 50 pc, n ~ 
10 2 cm -3 , and M ~ 10 5 -10 6 M ) and are typically larger 
and more massive than Bok Globules (D ~ 0.1-2 pc, n 
~ 10 3 -10 4 cm" 3 , and M - 1-100 M ? : IRathborne et al.l 
2006). Many have argued that stars can form anywhere 
in a GMC, but that the confining pressure and column 
density of an IRDC are requi red for the formation of clus- 
ters and massive stars (e.g. iMcKee fc Tad 12001 . The 
formation mechanism of an IRDC within a GMC is not 
well understood. Perhaps it is simply the gravitational 
collapse of a high-density perturbation in a GMC left on 
its own for many free-fall times. Alternatively, IRDCs 
may be the result of triggered collapse in a high-pressure 
environment, such as a cloud-cloud collision, or the dy- 
namic force from a bubble, such as an H II region, wind, 
supernova, or superbubble. Converging flows are another 
possible trigger for the collapse of an IRDC from part 
of a GMC. An IRDC, on ce formed, has the p roperties 
of cluster-forming clumps (|Lada fc L ada 2003) with the 
exception that IRDCs are much colder (10 - 20 K). Do 
these IRDC clumps fragment into cores which evolve in 
isolation or do the resulting cores undergo competitive 
accretion? This is a central quest i on in massive star 
forma tion (jZinnecker fc Yorkd 120071 IMcKee fc Ostrikerl 
l200l . 

This paper explores a sample of 17 IRDC clumps em- 
bedded within 8 IRDCs. We compare various mass trac- 
ers in IRDCs, distinguish active and quiescent clumps, 
and discuss a possible evolutionary sequence. This pa- 
per is organized as follows. In £j2j we describe the source 
selection and observations. In £]3J we present mass es- 
timation techniques, their results, and an analyis of the 
systematic errors involved. These masses are compared 
in §4.11 In §4.21 we compare the 8 fim dust extinction 
with the BGPS 1.1 mm dust emission. In £14.31 we ad- 
dress whether or not the IRDC clumps are bound, and 
in £14.41 we compare the molecular line tracers. In £14.51 
we discuss the radio continuum sources and stellar type 
limits for the UCH II regions found and in £|4.6l we discuss 
star formation activity in the clumps. In £j5j we discuss 
a possible evolutionary sequence for IRDC clumps. £j6]is 
a summary of our conclusions. Appendix A includes a 
discussion of each IRDC as well as images of each source 
with all the data used in the paper. 

2. SOURCE SELECTION AND DATA 

2.1. Source Selection 

In this paper, we investigate a sample of 17 IRDC 
clumps embedded within 8 IRDCs. Th ese were se- 
lected from the MSX dark cloud catalog (Simon et al. 
I2006a|) . and from the previous work of IRathborne et aL 



(2006), who selected the 38 highest contrast IRDCs from 
the dark cloud catalog with known kinematic distances. 
These 38 IRDCs, therefore , are among t he dar kest and 
densest in the catalog by iSimon et al.l (l2006aft. From 
this su bset, and based upon the work bv lChambers et all 
(2009), we selected mostly active clumps, with a few qui- 
escent clumps for comparison. Our sample consists of 
9 active, 4 intermediate, and 4 quiescent clumps. We 
mostly selected clumps with signs of active star forma- 
tion, motivated by our UCH II region survey, as this 
would maximize our detection rate. We included a few 
quiescent clumps in our sample for comparison, although 
no UCH II regions were expected in the quiescent clumps. 
Our sample of IRDC clumps is a subset of the darkest, 
most active IRDC clumps in the First Galactic Quad- 
rant (0 < 1 < 9 °) I n a sample of 106 IRDC clumps, 
IChambers et al.l ((2009) find that 65% are quiescent and 
35% are active clumps. We quantify the star-forming 
activity by four criteria: 1) an embedded 24 um p oint 
source, 2) "green fuzzies," ()Chambers et al.l |2009|) an 
extended 4.5 /im enhancement believed to be caused 
by shocked H2, thus indicative of ou tflows (also known 
as E xtended Green Objects, EGOs Cvganowsk i et all 
120081 ). 3) H 2 or CH 3 OH maser emission, and 4) Ultra- 
compact (UC) H II regions. We assign the designation 
"active" to clumps which have three or four of these signs 
of star formation, which means that each "active" clump 
has either an H II region or at least two outflow tracers 
and a 24 /um point source, and is actively forming stars. 
The designation "intermediate" is for clumps which ex- 
hibit one or two signs of active star formation. These 
clumps, therefore, may have only shock/outflow signa- 
tures or a 24 /im point source. We reserve the designa- 
tion "quiescent" for clumps with no signs of active star 
formation. 

Selection of the clump positions from the sample dis- 
cussed above required simply that there be an IR-dark 
feature, which can be near but not coincident with an IR- 
bright region, which has corresponding millimeter emis- 
sion. The exact clump positions were selected by eye 
based on the above criteria and the overlap with existing 
data. While the exact clump positions were not selected 
by stringent threshholds, visual inspection of the images 
verifies that the sample accurately represents the IRDCs. 
The clump p ositions and the c l osest millimeter clump, as 
identified bv IRathborne et al.1 £2006), are given in Table 

m 

2.2. Observations & Data Description 

For this study, we examine a sample of 17 IRDC clumps 
within 8 IRDCs. We utilize existing and new data to de- 
termine the physical properties and evolutionary stages 
of these clumps. A summary of the data used is found 
in Table d 

2.2.1. Mid- Infrared Data 

We utilize mid-IR data taken as part of the Galac- 
tic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE: iBenjamin et al.ll2003fl and 24 Mm data taken 
as part of the MIPSGAL survey (|Carev et al.|[2009h . Ex- 
tinction at 8 /im gives an independent mass estimate, the 
4.5 /im band determines the presence of "green fuzzies," 
and embedded stellar sources appear at 8 and 24 (im. We 
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TABLE 1 
IRDC Clump Positions 



IRDC Name 3 


Glimpse 


Bolocat b 


Closest MAMBO c 


R.A. 


DEC. 




Clump Name 


Clump Name 


Clump Name 


(J2000) 


(J2000) 


n /rcvnri f~~^ non oc i nn ai 

MoaDL. GU22.d5+UU.41 


GLM1 


2860 


MM1 


"i o on n a o 

lo:oU:24.2 


nn 1 n o o c\ 

-09:10:38.9 




GLM2 


2858 


TV /TTV \ O 

MM 3 


1 o on ot 

lc-:oU:o7.o 


nn in 1 

-1)9:12:54.1 


MbAJJU CjUzo.dU+UU.UU 


LxJjJVLJL 


3125 


MM 2 


1 0.Q/I.O1 Q 

lo:o4:21.o 


-Uo:lo:U(.o 




GLM2 


3132 


MM 7 


1 o o a nn 

18:34:26.6 


no "i T ni t 

-08:17:21.7 


MoADG GU24.ooH-UU.il 


GLM1 


3284 


MM1 


"i o o tr no 1 

lo:o5:Uo.l 


nT or m n 

-U7:oo:U1.9 






onon 
OZoZ 


A/TA/T/1 
1V11V14 




fl7.Q7. in 7 


MSXDC G024.60+00.08 


GLM1 


3382 


MM2 


18:35:35.8 


-07:18:22.2 




GLM2 


3383 


MM1 


18:35:39.9 


-07:18:46.0 


MSXDC G028.23-00.19 


GLM1 


3923 


MM1 


18:43:30.3 


-04:13:03.7 


MSXDC G028.37+00.07 


GLM1 


3946 


MM14 


18:42:53.4 


-04:02:23.8 




GLM2 


3939 


MM4 


18:42:50.6 


-04:03:16.9 




GLM3 


3942 


MM6 


18:42:48.9 


-04:02:05.4 




GLM4 


3955 


MM1 


18:42:52.5 


-03:59:48.1 


MSXDC G028.53-00.25 


GLM1 


3998 


MM2 


18:44:16.1 


-04:00:09.9 


MSXDC G034.43+00.24 


GLM1 


5373 


MM1 


18:53:18.0 


+01:25:24.9 




GLM2 


5373 


MM2 


18:53:18.4 


+01:24:51.3 




GLM3 


5385 


MM3 


18:53:19.6 


+01:28:24.6 



a from the MSX dark cloud catalog by Simon ct al. ( 2006a]). For the remainder of the paper we drop the MSXDC 
prefix. 

BGPS 1.1 mm clump source catalog, Bolocat, number { Rosolowsky ct aT1 l20107) 
c from lRathborne et all .20061) 



TABLE 2 

Summary of Observations 



Data 


A 


Beam FWHM 


la Sensitivity 


GLIMPSE 


3.6 fim 


1.7" 


0.3 MJy/Sr a 


GLIMPSE 


4.5 /an 


1.7" 


0.3 MJy/Sr 


GLIMPSE 


5.8 fim 


1.9" 


0.7 MJy/Sr 


GLIMPSE 


8.0 )im 


1.9" 


0.6 MJy/Sr 


MIPSGAL 


24 /an 


6" 


0.67 mjy 


HHT 


1.12 mm b 


28" 


0.04 K 


HHT 


1.07 mm b 


27" 


0.04 K 


GRS 


2.86 mm b 


46" 


0.27 K 


VLA 


3.6 cm 


2.8 x 2.4" 


0.05 mjy/beam 


IRAM 


1.2 mm 


11" 


10 mjy/beam 


BPGS 


1.1 mm 


33" 


30 mjy/beam 



a GLIMPSE raw la surface brightness sensitivities from 
IReach et ail 1 12009) 

b HCO+ J = 3-2, N 2 H + J = 3-2, and 13 CO J = 1-0 respec- 
tively. For HCO+ and N 2 H+, 0.04 K is the lcr sensitivity per 
1.1 km s _1 channel and for 13 CO, 0.27 K is the la sensitivity 
per 0.21 km s _1 channel 



utilize the "green fuzzy" (ext ended 4.5 ixm emission , in- 
dicative of shocks) catalog bv lChambers et all ()2009l ) . as 
well as their 24 /im point source identification and H2O 
and CH3OH maser survey. 

2.2.2. Millimeter Continuum Dust Emission 

The millimeter continuum data traces the cold dust 
within the IRDCs and gives cloud and clump mass esti- 
mates, as well as the locations of compact clumps within 
the IRDCs. We utilize 1.2 mm data obtained at the In- 
stitut de Radioastronomie Mill imetrique (IRAM) 30 m 
telescope with MAMBO II from lRathborne et all (|2006l) . 
The FWHM angular resolution of each element in the ar- 
ray is 11" , with a separation between each element of 20" . 
The final (1 a) r.m.s. noise level in these maps is ~ 10 
mjy/beam. 

We present additional millimeter data taken with Bolo- 
cam (|Glenn et al.l 120031) on the Caltech Submillimeter 



Observatory (CSO) 10 m diameter telescope at 1.1 mm 
as part of the Bolocam Galactic Plane Survey, BGPS vl.O 
PI (jAguirre et al.ll2010HRosolowskv et al.ll2010f ). Bolocam 
is a 144-element (of which 115 are working) bolometer ar- 
ray arranged on a uniform hexagonal grid. Each bolome- 
ter has an effective FWHM of 33" and the array field 
of view is 7 '.5. The data were taken with a 45 GHz 
bandwidth filter centered at 268 GHz, which excludes 
the bright 230 GHz J=2-l CO line. The data were re- 
duced using a custom pipelin e created for the BGP S us- 
ing the methods described in lAguirre et al.l ([2010). The 
dominant source of noise at these wavelengths, the at- 
mospheric component, is removed using Principle Com- 
ponent Analysis (PCA). This cleaning process limits the 
spatial scale to which we are sensitive to between about 
30" to 400". The final (1 a) r.m.s. noise level in these 

1 http://irsa.ipac.caltech.edu/data/BOLOCAM_GPS/ 
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maps is ~ 30 mJy/beam. 



2.2.3. 



CO J =1-0 Molecular Gas Tracer 



As a tracer of column density, radial velocities, and for 
spectral line comparison, we include 13 CO J=l-0 data 
taken as part of the Boston University-Five College Ra- 
dio Astronom y Obsevatory Galact ic Ring Survey (BU- 
FCRAO GRS; I Jackson et all [20061) . These data was ob- 
tained from 1998-2005 using the single sideband focal 
plane array receiver, SEcond QUabbin Optical Imaging 
Array (SEQUOIA) on the FCRAO 14 meter. At a cen- 
tral frequency of 110.2 GHz, the beam FWHM is 46" 
and the spectral resolution is 0.21 km s _1 . We correct 
the 13 CO antenna temperatures for the main beam ef- 
ficiency of 0.48 through the expression T mb = T\jr\ mh . 
The 1 a sensitivity in the GRS data is about 0.27 K. 

2.2.4. VLA Radio Continuum Data 

In order to unambiguously determine which massive 
clumps have entered the main sequence, we used the Very 
Large Array (VLA) in New Mexico to look for thermal 
bremsstrahlung from newly formed H II regions. We ob- 
served a total of twelve positions over eight IRDCs in 
the X band in the C configuration in the spring 2008. 
We observed single channel continuum in the X band, 
centered at 3.6 cm with a bandwidth of 50 MHz. We 
performed flux calibrations on 3C286 and phase calibra- 
tions on 1832-105 and 1851+005. The data were reduced 
using standard continuum interferometric techniques in 
the AlPS software package. The primary beam is 5.4', 
and the highest resolution element is about 2.8 x 2.4". 
This translates to 0.06 (0.10) pc resolution for the nearest 
(farthest) clump. We had a total of three six-hour ses- 
sions. We observed each source three times per session, 
changing the order each night to obtain good coverage 
of the u-v plane. We integrated for an hour total on 
each source achieving a 1 a r.m.s. flux density of ~ 0.05 
mJy/beam. 

2.2.5. HCO + and N 2 H + Dense Molecular Gas Tracers 

To determine the radial velocities of the clumps and 
their molecular structure, we probe the dense molecu- 
lar gas using the J=3-2 transitions of HCO + and N2H+ 
(rest frequencies of 267.5576190 and 279.5117010 GHz 
respectively) . Using the 1 mm dual polarization ALMA- 
prototype Sideband Separating Receiver on the Hein- 
rich Hertz Telescope (HHT), we observed HCO + and 
N2H + simultaneously with 1.1 km s _1 spectral resolu- 
tion and 28" and 27" spatial resolution (O fwhm)- We 
performed two minute position-switched observations on 
the seventeen IRDC clumps, which gave us an r.m.s. 
noise level of about .04 K. We On-the-Fly mapped (see 
iMangum et al. 1 120071) the five most compelling clouds in 
both molecular lines simultaneously. We mapped the 
clouds, scanning twice in RA and twice in DEC with 
a row spacing of 10" at a scan speed of 5" per second. 
The separate maps were offset by about 4" . We pointed, 
focused, and calibrated on Jupiter every 2 or 3 hours 
when it was available, otherwise, we pointed and focused 
on G34.3+0.15. We used S140 as a check of the line 
position and intensity about every two hours, and also 
to determine the Vertical / Horizontal Polarization ratio 



(Vpol/Hpol) and the typical sideband rejection measure- 
ment. Typical sideband rejection between the USB and 
LSB for each polarization was between 13 and 15 dB. We 
correct for the main beam efficiency of 0.75 and the ratio 
Vpol/Hpol = 1.1 in our data reduction. 

The single spectra were reduced using standard pro- 
cedures in the Gildas CLASSQ software package. The 
baselines were removed, the horizontal polarization was 
multiplied by the Vpol/Hpol scaling factor and then co- 
added. We then corrected for the main beam efficiency, 
giving us T m b. The On-the-Fly maps were reduced using 
some custom and shared software in the Gildas CLASS 
package. The same procedure as above was followed 
for each individual spectrum. The spectra were then 
interpolated onto a regular grid using a Gaussian ta- 
pered Bessel functio n (Reiter, M. private communication; 
IMangum et"aTll2007t ). 

3. MASS ESTIMATES 

Mass is one of the most fundamental physical proper- 
ties necessary to understand the nature of a region. How- 
ever, there is no perfect tracer of mass. Molecular line 
mass tracers, such as 13 CO may suffer from depletion, op- 
tical depth effects or varying excitation conditions. Virial 
masses are only reliable in the idealized case of a gravita- 
tionally bound object in virial equilibrium, where the line 
is broadened simply by virialized motion. An extinction 
mass is limited by a priori knowledge of the foreground 
and background radiation field, and the absorption prop- 
erties of the obscuring dust. Dust continuum massses are 
thought to be a promising tracer of cold, dense environ- 
ments, but are also sensitive to varying dust temperature, 
emissivity, and grain properties. The distance estimate 
plays a significant role in the mass uncertainty via the 
assumption of a Galactic rotation curve, and the ambigu- 
ity associated with the velocity scatter of objects at any 
given distance. Each of these tracers is also affected by 
the statistical uncertainties in each measured quantity. 
For this study, we find that the systematic uncertainties 
far outweigh the statistical uncertainties and include only 
estimates of the systematic uncertainty in §3.61 

A minimization and characterization of systematics is 
vital for a fair comparison of IRDC mass tracers. For this 
reason, we have convolved the extinction and MAMBO 
1.2 mm maps to the BGPS resolution using a Gaus- 
sian kernel with a width given by a conv = J a% gps - a%. g 

where a org is the original width of the beam that is being 
convolved. Each mass is then retrieved within a circular 
aperture with the equivalent area of the BGPS Gaussian 
beam, 

O fwhm 

r *ff = IT, 



= 19.81". 



(1) 



This allows us to compare these mass estimates indepen- 
dent of beam filling factors and beam or source sizes. 
Where possible, we also calculate the total cloud mass, 
the mass within the elliptical aperture used to define the 
extent of the IRDC for extinction mapping (see Figure[I|. 
The line widths used to determine the virial masses were 
achieved from single pointings of HCO + on the HHT, 
which has an effective beam FWHM of about 28" at these 

1 http://www.iram.fr/IRAMFR/GILDAS 
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TABLE 3 
Properties of the IRDC clumps 















J co 


HCO + 


N2H + 


Integrated Intensity (K km s 1 ) 


Clump 


Vlsr 


A V 


Distance 


fforc 




Peak T MB 


Peak T MB 


Peak T MB 


1 CO 


HCO + 




(km s 1 ) 


(km s -1 ) 


(kpc) 




(K) 


(K) 


(K) 








G022. 35+00. 41 
























GLM1 


52.4 


3.4 


3.6 


0.14 


1.3 


5.4 


1.9 


1.2 


8.6 


8.0 


5. 1 


GLM2 


83.9 


3.4 


4.8 


0.23 


0.7 


8.2 


0.8 


0.3 


18.0 


2.5 


0.8 


G023.60+00.00 
























GLM1 


53.3 


6.7 


3.6 


0.14 


0.7 


4.3 


1.1 


1.3 


12.0 


7.1 


6.0 


GLM2 


53.4 


4.5 


3.6 


0.14 


0.8 


4.7 


0.6 


0.6 


12.1 


2.4 


2.4 


G024. 33+00. 11 
























GLM1 


113.4 


6.7 


5.9 


0.33 


0.7 


6.8 


1.6 


3.4 


30.0 


12.1 


20.0 


GLM2 


114.1 


5.6 


5.9 


0.33 


0.7 


5.2 


0.5 


0.3 


26.1 


2.8 


1.8 


G024. 60+00. 08 
























GLM1 


114.7 


3.1 


6.0 


0.33 


1.3 


4.6 


0.6 


0.5 


7.7 


2.2 


2.2 


GLM2 


53.0 


5.6 


3.5 


0.14 


1.2 


•5.6 


1 . 5 


1.0 


14.1 


8.7 


5.5 


G028.23-00.19 
























GLM1 


79.7 


1.5 


4.6 


0.23 


1.5 


o.8 


0.3 


0.3 


28.4 


1.7 


1.5 


G028.37+00.07 
























GLM1 


78.8 


5.6 


4.5 


0.23 


2.6 


4.7 


1.1 


1.2 


21.0 


5.9 


4.1 


GLM2 


77.9 


4.5 


4.5 


0.23 


2.5 


6.3 


2.1 


1.6 


29.0 


9.9 


6.8 


GLM3 


78.9 


2.2 


4.5 


0.23 


1.8 


4.9 


1.2 


1.3 


24.7 


4.9 


5.2 


GLM4 


76.7 


6.7 


4.4 


0.23 


1.4 


7.0 


1.4 


1.8 


25.9 


7.4 


12.3 


G028.53-00.25 
























GLM1 


86.2 


5.6 


4.9 


0.26 


0.4 


6.6 


0.4 


0.2 


22.3 


1.7 


1.1 


G034.43+00.24 
























GLM1 


58.0 


6.7 


3.6 


0.18 


0.7 


6.7 


5.1 


6.1 


23.9 


37.4 


37.0 


GLM2 


57.5 


6.7 


3.6 


0.18 


0.4 


8.5 


6.4 


6.2 


30.0 


46.0 


34.6 


GLM3 


59.5 


5.6 


3.7 


0.18 


0.9 


7.2 


2.9 


2.5 


18.4 


18.0 


11.7 



a The fraction of foreground emission along the line of sight, see i|3.2l 

frequencies, similar to the BGPS beam. In most cases 
the radius used to determine the virial mass (the "virial 
radius" ) is very similar in size to the circular aperture ra- 
dius, so these are fair comparisons in terms of source size. 
The 13 CO masses were retrieved from single GRS point- 
ings, with an effective FWHM of about 46" . The mass in 
this case is estimated over a greater area than the other 
techniques, so it will include a wider region and more 
emission, which we consider when making mass compar- 
isons. We chose not to degrade the resolution of the 
other maps past the 33" for the BGPS beam to match 
the 13 CO beam as the 13 CO is considered a secondary 
mass estimator, and is included only for comparison. 

We explain below the techniques for achieving the dif- 
ferent mass estimates and the estimated systematic un- 
certainty in Section 13.61 We then compare these mass 
estimates and discuss their differences and similarities. 

3.1. Distance Estimates 

CO is the second most abundant molecule in inter- 
stellar space. It is found nearly everywhere, includ- 
ing IRDCs. In these dense, cold environments, how- 
ever, even the less abundant 13 CO can bec ome optically 
thick and can deplete onto dust grains fe.g. lTafalla et al.l 
l200l . HCO+ and N 2 H+ are less abundant than 13 CO 
and therefore become optically thick at higher columns 
of H2. The critical density is also much higher, 3.5 x 
10 6 cm -3 and 3.0 x 10 6 cm" 3 for the J=3-2 transitions 
of HCO + and N2H + respectively, compared with 1.9 x 
10 3 cm -3 for 13 CO J=I-0. Therefore, there is typically 
only one HCO + or N2H+ source along a given line-of- 
sight, which can then be confidently associated with the 
dense IRDC, as opposed to CO which is often seen at 
many velocities along a given line-of-sight through the 



Galactic plane. For this reason, we use the HCO + veloc- 
ity to determine the kinematic distance to each IRDC. 
A kinematic distance had already been determined to 
eac h of these ob j ects us ing 13 CO morphology matching 
by iSimon et al.l (|200 6bj. however, we find that in the 
confused inner Galaxy, morphology matching can mis- 
assign distances. Two clouds of our sample of eight 
(G022.35+00.41 and G024.60+00.08) are spatially adja- 
cent clumps along the line-of-sight at very different dis- 
tances, as opposed to a contiguous object as found by 
ISimon et al.f(|2006bf ) . 

Distances were estimated kinematically using the 
HCO + clump v elocities and the rotation curve of 
iReid et all (|2009f l . determined with trigonometric paral- 
laxes and proper motions of masers using the VLBA and 
Japanese VERA project. In the First Galactic quad- 
rant, a radial velocity can correspond to two possible 
distances, a problem known as the "Kinematic Distance 
Ambiguity" . Since IRDCs are perceived as dark extinc- 
tion features against the Galactic background, it is as- 
sumed that all IRDCs are located at the near kinematic 
distance. If we assume an average velocity scatter at 
any given distance of about 5% (about 5-10 km s _1 ), 
this translates to an approximate distance uncertainty 
of 20%. The distances, line-widths, and observed prop- 
erties of the clumps are listed in Table |3] 

3.2. Extinction Masses 

The mid-IR extinction by IRDCs provides a power- 
ful way to obtain an independent measure of their col- 
umn density. We restrict our analysis of extinction to 
the 8 fim band of GLIMPSE, as this band is dominated 
by emission from the diffuse ISM, rather than stellar 
sources. We apply the extin ction mapping method put 
forth bv lButler fc Tanl (|2009h . In brief, this method esti- 
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Fig. 1. — G023. 60+00. 00 Top: 8 /im GLIMPSE image with ellipse approximating the extent of the IRDC. This image represents the 
values of l v i :0 b s in Equation [5] Middle: Circularly median filtered 8 (im image with the points in the ellipse determined by interpolation 
from outside so they are not skewed by the IRDC. This is a proxy for the diffuse Galactic background, I„q i, s in Equation [4] Bottom: Final 
surface density map in units of g cm -2 



mates the Galactic background and foreground emission, 
finds the optical depth required to produce the observed 
extinction feature, assumes a dust opacity, and finally 
extracts the gas mass surface density. Our application of 
this method is described below. 

To accurately map the extinction, we need to estimate 
the intensity of radiation behind the cloud of interest, 
I u o, and the radiation in front of the cloud, \ v \. As- 
suming no emission from the cloud itself at these wave- 
lengths, we have 

hi = Loe~ Tv (2) 

where the dust opacity, and £ is the 

surface mass density. We adopt a value of Ks^m = 11-7 
cm 2 g" 1 , following the analysis of iButler fc Tanl (|2009f ) 
(they adopt a value of 7.5 cm 2 g _1 as they use a gas- 
to-dust ratio of 156, rather than the value of 100 used 
in this paper). This value is closest to the model of 
lOssenkopf fc Henninsj (|1994f ) of thin ice mantles that 
have undergone coagulation for 10 5 years at a density 
of nn 2 ~ 10 6 cm -3 , which is a reasonable model for the 
cold, dense environment of an IRDC. This expression as- 
sumes that the dust opacity remains constant along the 
line-of-sight, which may not be true. However, most of 
the extinction arises from the coldest, densest parts of 
the cloud. In the less dense regions, the opacity may 
vary by up to a factor of two, but we assume that the 
majority of the extinction is caused by the coldest, dens- 



est regions of the cloud, which are well-characterized by 
the adopted value of /cs^m- 

Before proceeding with the background determination, 
we must first estimate the fraction of emission from 
dust in the foreground. We assume that the Galactic 
distribution of hot dust (in the Galactic midplane) is 

iven by the Galactic surf ace density of OB associations 

McKee fc Williamsl[l997l ) 



Sob oc e~ R / H « 



(3) 



where R is the Galactocentric radius (in kpc) and 
= 3.5 kpc is the Galactic radial scale lengt h We then 
integ rate the column of dust from the sun (jReid et al.l 
2009, at Ro = 8.4 kpc from the Galactic center) to the 
cloud, and the column of dust along the same line of 
sight, from the Sun out to a Galactocentric radius of 16 
kpc (beyond which there is a negligible contribution to 
the dust emission) . The ratio of the column to the cloud 
over the column out to 16 kpc is called the 'Foreground 
Intensity Ratio', ff orc ; the fraction of emission along the 
line-of-sight that is from the foreground. Thus, the true 
radiation behind the cloud, I„0j determined from that 
observed, I„o,o6s, is given by 

IvO = (1 — ffore)lu0,obs (4) 

and the true radiation in front of the cloud, I„i, deter- 
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mined from that observed, b s is given by 



I, 



Iu 



Lvl,obs f fore^uO^obs • (5) 

iButler fc~T an (200S|) note that ffore is uncertain because 
of the inevitable small-scale variations in the Galactic 
background and the fact that the regions surrounding 
the IRDC are likely to be embedded in the same GMC 
that hosts the IRDC. This leads to a higher extinction of 
the integrated Galactic background, creating a tendency 
to underestimate ff ore and thus E. 

An upper limit for ff ore is provided by the minimum 
flux at 8 /im divided by the average background. This 
value would be ff ore if the cloud were totally opaque at 
8 /im. This upper limit is about 2 times higher in most 
cases, except in G028. 37+00. 07, where it is about the 
same. The remainder of the IRDCs may or may not be 
optically thick at 8 /im, however, the good correlation 
(see Figure [2]) with BGPS 1.1 mm masses (which is al- 
most certainly optically thin) is an indication that many 
of the IRDCs are still optically thin at 8 /zm. 

For the purpose of estimating the Galactic background, 
we employ the more a ccurate small- s cale m edian filter 
(SMF) approach from IButler fc Tanl (pOOl ). This ap- 
proach is designed to capture the small-scale variations 
in the background, without altering the estimate by the 
darkness of the IRDC itself. Thus, we first exclude an 
ellipse approximating the size and shape of the IRDC 
from the filtering process (Figure [T] top). To smooth the 
data, we first remove the high-end tail of the distribution 
(all points greater than twice the mode) which is mostly 
due to stars. We then perform a circular spatial median 
filter at each point, i.e. we compute the median of all the 
data in a circle around a given point, and that becomes 
the smoothed data value at that point. A reasonable size 
for the filter was empirically determined to be a radius 
of one arcminute. This smoothed map is our estimate 
of the diffuse Galactic background. We then estimate 
the 'background' inside the ellipse by taking the aver- 
age of the smoothed data values from outside the ellipse, 
weighted by the inverse square of their distance from the 
point. This is depicted in the middle frame of Figure [TJ 
This image is now our estimate of Ij/o, obs, and our orig- 
inal image is bs- Using Equation [5] we then create 
a surface density map, shown in Figure [T] in the bottom 
panel. 



1 



-hi 



fl.obs 



fforcI^O.obs 



(1 — ffore)li/0 5 obs 



(G) 



There is an additional correction that must be 
taken into account for quantities derived from ex- 
tended sources in the IRAC data. Scattering from 
extended sources within the IRAC array systemtically 
increases the measured surface brightness. For a 
uniformly distributed extended source, this scattering 
increases the measured surface brightness by about ~ 
30% at 8 lira ( S. Carey, private communication and 



by Ti = Li.obs + s I„o,o& s - Plugging these values into 
Equation |6l simplifying, and writing the surface density 
in terms of our measured quantities, our final expression 
for the surface density is: 



£ = -—In 



(S + l)I„l,obs — (S + f)fforcliy0,obs 
(1 — fforc)lj/0,obs 



(7) 



Using distances d e termin ed from the HCO + velocities 
and the IReid et all (|2009D rotation curve, we are then 
able to sum the surface density over the BGPS beam- 
sized apertures to obtain the mass (see Table 0] and 
Figure in the clump. The extinction masses for the 
clumps range from 60 to 620 M Q and from 1200 to 26,020 
M Q for the whole IRDC. 

3.3. Dust Continuum Masses 

Obtaining mass estimates of IRDCs with molecular 
lines is often problematic due to high column densities 
and molecular freezeout onto dust grains. The cold dust 
emission, however, as observed in the millimeter, pro- 
vides more reliable mass estimates as it is optically thin 
and does not deplete. We use both the 1.2 mm and 
1.1 mm dust continuum emission from MAMBO and the 
BGPS, to achieve isothermal mass estimates using the 
expression 

M = S " D * s (8) 

where S„ is the source flux density, D is the distance, 
k v is the dust opacity and B y (T) is the Planck func- 
tion at dust temperatu re, T. We adopt a va lue of Ki.i mm 
= 0.0114 cm 2 g~ x from lEnoch et al.l (|2006l ). in which we 
have assumed a gas -to-dust ratio of 100. This op acity is 
consistent with the lOssenkopf fc Hcnning ( 1994) model 
used for the 8 /zm opacity in the extinction masses. The 
mass equation reduces to 



M 1Amm = 14.32 (e 13 °/ T - 1) ( A- ) ( JL- ) M 



1 Jy / \ 1 kpc 



(9) 



M M = 20.82 (e^-l)(A) j^j M . 

(10) 

We assume a dust temperature T = 15 K. We can also 
find the column density of H 2 , as given by 



N(H 2 ) 



S v 



VL B KvB u (T)n H2 m H 



(11) 



N(H 2 ) 1 



= 2.20 x 10 22 (e 13 °/ T - 1)S„ cm" 2 (12) 



http://ssc.spitzer.caltech.edu/irac/iracinstrumenthandbook/33^^ 2 ) 
This means that the measured background (l^o, obs) is 
systematically brighter than the true background, 
which we will call To for simplicity. Where s = 0.3, 
the scattering correction, the true background is given 
by To = (1 + s)I^o, obs- The measured foreground 
(Ij/i, obs) is also contaminated by scattered light from 
the background, so the true foreground (Ti) is given 



3.20 x 10 22 (e 12 °/ T 



1)SV cm" 2 (13) 



where £Ib is the beam size, /i# 2 is the mean molecu- 
lar weight for which w e adopt a value of /i# 2 = 2.8 
(|Kauffmann et al.l [2008) . The dust continuum masses 
and column densities are listed in Table @] The total 
BGPS 1.1 mm cloud masses range from 1200 to 13,000 
M Q while the clump masses range from 80 to 1100 M Q . 
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TABLE 4 

Mass Comparison 





Aperture 


Column Density in 33' 


BGPS Beam 


1J CO 


Mass inside 33" BGPS Beam 


"CO 






Clump a 


Size b 


Extinction 


BGPS 


MAMBO 


Column 


Extinction 


BGPS 


MAMBO 


Mass 


Virial Mass 






(PC) 




(10 22 cm- 


2 ) 




(M ) 


(M ) 


CM©) 


(M ) 


(M ) 


(PC) 


CjU22.do+UU.41 
























GLM1 


0.7 


2.3 


1.1 


2.9 


0.4 


190 


95 


250 


70 


340 


0.2 


GLM2 


0.9 


1.3 


1.2 


1.4 


0.9 


200 


180 


210 


240 


2300 


1.3 


CjUzd.oUH-UU.UU 


2.9 x 1.8 










1200 


1200 


1800 








GLM1 


0.7 


1.2 


1.7 


3.0 


0.6 


100 


1 10 


250 


100 


1800 


0.3 


GLM2 


0.7 


1.4 


1.0 


1.7 


0.6 


120 


81 


140 


90 


660 


0.2 


CjU24.Jd-l-UU.ll 


6.0 x 8.0 










12010 


13000 


14000 








LrL/JVLl 


1. 1 


1.4 


4.8 


8.6 


1 . 5 


310 


1200 


1950 


610 


2500 


0.1 


GLM2 


1.1 


1.3 


1. i 


1.9 


1.3 


300 


310 


440 


530 


3100 


0.7 


CjU24.oU+UU.Uo 
























GLM1 


1.2 


2.4 


1.3 


2.3 


0.4 


560 


310 


550 


160 


850 


0.5 


GLM2 


0.7 


2.1 


1.6 


3.2 


0.7 


170 


130 


260 


100 


1000 


0.2 


G028. 23-00. 19 


5.2 x 3.3 










5520 


3700 


1900 








O T T\ A 1 


0.9 


2.8 


1.1 


2.0 


1 . 1 


380 


190 


270 


350 


910 


0.3 


G028.37+00.07 


5.8 x 4.5 










26020 


13000 


18000 








GLM1 


0.9 


4.7 


1.6 


3.0 


1.0 


620 


220 


400 


250 


2100 


0.5 


GLM2 


0.9 


4.6 


2.2 


3.9 


1.4 


610 


290 


520 


350 


330 


0.1 


GLM3 


0.9 


3.3 


1.4 


2.6 


1.2 


430 


190 


340 


290 


370 


0.5 


GLM4 


0.9 


2.6 


5.3 


10.3 


1.3 


350 


700 


1400 


290 


4400 


0.7 


G028.53-00.25 


7.3 x 4.6 










4800 


7300 


11000 








GLM1 


0.9 


0.7 


1.6 


3.5 


1.1 


110 


250 


540 


310 


3500 


0.8 


G034.43+00.24 


6.9 x 2.3 










5550 


7500 


14000 








GLM1 


0.7 


1.3 


10.1 


17.8 


1.2 


110 


850 


1500 


180 


3000 


0.4 


GLM2 


0.7 


0.7 


9.3 


15.9 


1.5 


60 


790 


1300 


230 


2500 


0,1 


GLM3 


0.7 


1.7 


2.2 


3.9 


0.9 


140 


190 


330 


150 


1700 


0.4 



a The clump names are from the MSX dark cloud catalog ]Simon et al.ll2006at) , and have had the "MSXDC" prefix removed 

k For the clumps, the aperture size is the physical size of the 40" diameter circular BGPS aperture. For the cloud, it is the physical major and minor 
axes of the ellipse approximating the extent of the IRDC. 

3.4. Virial Masses 

The molecular line width of these objects provides 
an independent measure of mass. The virial mass is 
based on the assumption that the object is gravitation- 
ally bound and in virial equilibrium. Therefore, the 
virial mass estimate, in combination with other mass 
estimates, may allow discussion of whether or not the 
clump is bound and virialized. For a density profile of 
p(r) = r _Q , the virial mass is given by 



2a 



3 — a 



Ra 2 
G 



(14) 



where a is the line-of-sight velocity dispersion, R is the 
virial radius, and G is the gravitational constant. We 
adopt a spectral index of a = 1.8, which was found to be 
th e mean in a sample o f 31 massive star-forming regions 
bv lMueller et al.l ()2002l ) with a standard deviation of 0.4. 
However, this spectral index was derived from a sample of 
"active" clumps, so this is an added source of uncertainty 
for quiescent clumps. We do not include a correction for 
ellipticity, as these clumps are compact enough to be 
well- approximated by spheres. This expression reduces 
to 



M vir = 147 



where R is the virial radius and Avf wn m is the line width 
at half maximum. In order to determine the best es- 
timate for the virial radius, a two-dimensional circular 
Gaussian was fit to each BGPS clump. The Gaussian fit 



was then deconvolved from the BGPS beam. The effec- 
tive radius (see Eq. [I} of the deconvolved Gaussian fit is 
the virial radius, given in Table[4] In this study, we assign 
a virial 'radius' to a Gaussian beam which is likely con- 
volved with a power-law distribution of cores within the 
beam. The meaning of the virial 'radius' in this case is 
tricky as the deconvolved source size can vary by up to a 
factor of two depending on the slope o f the power law dis- 
tribution of sources within the beam (|Shirlev et al.|[2003l 
Figure 16). However, this method is consistent between 
the clumps in our sample, and is a reasonable proxy for 
the size scale of the diffuse emission. 

We use the FWHM of the HCO+ line for the virial 
masses presented here. Using the N2H+ line width in- 
stead does not change the results qualitatively, as the two 
tracers provide consistent line widths across different ac- 
tivity levels. N2H+ has been shown to be well correlated 
with submm dust column in starless cores (|Tafalla et al.l 
120021 : IWvrowski et al.l [2000). meaning that these virial 
masses are most likely based on the same material which 
is producing the dust continuum. The virial masses of 
the clumps were found to range from 330 to 4400 M Q . 

3.5. 13 CO Masses 

13 CO is commonly used to estimate the column density 
and mass of molecular clouds in th e Galaxy. However , 
due to depletion onto dust grains dTafalla et al.l l20f)l 
and optical depth effects in cold, dense environments, 
13 CO is not the most reliable tracer of column in IRDCs. 
While there is good morphological agreement between 
the 8 (an extincti on and the 1 3 CO toward the dense 
clumps in IRDCs, IDu fc Yand (|200cl show that 13 CO 
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Fig. 2. — A comparison of mass estimates for the clumps (see Table[4]l. The line drawn has a slope of 1, representing equality. The red 
crosses are "active" clumps, the green are "intermediate," and the blue are "quiescent." Note that the four highest BGPS 1.1 mm mass 
objects (where there is significant discrepancy) are all associated with UCH II Regions. 



is optically thick toward many IRDCs, so we caution the 
reader that these tracers are probabl y not probing to the 
same depth in dense IRDC clumps. IDu fc Yand (poM ) 
also demonstrate that typcial excitation temperatures of 
13 CO in IRDCs are of the order of - 10 K. This indi- 
cates that the 13 CO is at the very least tracing to some 
depth in the cold, dense region of the clump, however the 
13 CO is also tracing the surrounding GMC and IRDC en- 
velope. We include column density and mass estimates 
from 13 CO as a secondary mass tracer, and as a compar- 
ison to other tracers. In the optically thin, thermalized 
limit, the 13 CO column density is given by 



expression then reduces to 

N(H 2 ) = 1.45 x 10 17 - 



jT mb dv 



(17) 



\ _ e -5.29/T 

This column density can then be converted to a mass via 

M = N{H 2 )^ H2 m H A (18) 

where fj>H 2 is the mean molecular weight of H2, m# is 
the mass of hydrogen, and A is the physical area over 
whi ch the mass is summed. We adopt a value of fiH 2 = 
2.8 ([Kauffrnann et aL l[2008). This expression reduces to 



-hv/kT c: 



T m bdv (16) 



where v is the frequency of the 13 CO J=l-0 transition, 
Am is the Einstein A coefficient of 13 CO from state J=l 



10 

to J=0, T ex is the excitation temperature, Bj is the 
rotation constant, and Xi3 C p is the fraction of 13 CO 
to H 2 . We adopt a value of 12 CO / 13 CO of 58 from 
ILucas fc Liszl (Il998h . and a value of 12 CO / H 2 of 1CT 4 , 
a value of 55.101038 GHz for Bj, and standard NIST val- 
ues for all constants and spectral transition values. This 



A 



D 



M = 5.00 x 10- 25 N(H 2 ) [ T ^—; ) ( — ] M 



1 kpc 



(19) 

where A is the area in arcsec 2 , over which the mass is 
summed, and D is the distance in kpc. The 13 CO clump 
masses range from 70 to 610 M Q . 

3.6. Mass Error Analysis 

For this study, we find that the systematic uncertain- 
ties far outweigh the statistical uncertainties. We present 
below our estimated systematic uncertainties and the 
method used to derive them. 
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We assign each variable in the mass equation (Eqns. 
El E3 Eil HS1 an d HU) a fractional uncertainty, perform a 
Monte-Carlo simulation randomly drawing the variables 
from a Gaussian distribution, and report the 68% confi- 
dence interval over which the mass varies. Since we have 
reason to believe that each of the variables is centered at 
our adopted value, we assign each variable a randomly 
populated Gaussian distribution, centered at the adopted 
value and with a width such that 99.7% of the points 
are within the lower and upper bounds of the fractional 
uncertainty. For example, we assign an uncertainty to 
the distance of 20%. Therefore, for each run, we ran- 
domly sample the distance value from a Gaussian which 
is centered at the adopted value and has a width such 
that 99.7% of the points lie between 0.8 and 1.2 times 
the adopted distance value. For each mass estimate, we 
ran at least 10 7 points in the simulation. We note that 
several of the mass probability distributions below are 
asymmetric. To be consistent with the literature, we 
present the mass estimates at the adopted variable val- 
ues (which is equal to the mean of the mass probability 
distribution), but note that this is not the same as the 
median or most probable value of the distribution. The 
68% confidence intervals were drawn from the median of 
the distribution outwards. We divide the high and low 
mass ends of the confidence interval by the mass at the 
peak of the distribution to achieve the asymmetric 68% 
confidence intervals. 

The dominant sources of uncertainty in the extinction 
mass estimate are the background and foreground esti- 
mates, the scattering coefficient, as well as the distance, 
opacity, and ff or c- We find that typical fluctuations in the 
8 /im emission near IRDCs, on the scale of the IRDCs, 
is around 40%. This background fluctuation is the best 
measure we have for how much the background behind 
the IRDC will vary from what we measure, so we assign 



a 40% uncertainty to the background estimate. In cal- 
culating the foreground uncertainty, we assign the same 
40 % uncertainty to the measured background and an 
additional 50% uncertainty to the value of ff or c- We in- 
clude a 15% uncertainty in the value of the scattering 
coefficient as the color variation is about 5 % and not 
being able to characterize the spatial dependency of the 
scattering has about a 10-15% effect (S. Carey, private 
communication). We also include a 20% uncertainty in 
the distance determination (as discussed in Section [Q]) 
and factor of two un certainty in the value of the o pacity 
as recommended by lOssenkopf fc HenningI (fi"994f) . The 
68% confidence intervals for the extinction mass are from 
80% to 2.2 times the quoted mass. 

For the dust continuum (BGPS 1.1 mm and MAMBO 
1.2 mm) mass estimates, the dominant uncertainties 
are the temperature, the flux density, the opacity, and 
the distance. For both the 1.1 mm and 1.2 mm 
data, we assign an uncertainty of 50% to the temper- 
ature ( for most clumps the temperature is between 10- 
20 K; iDunham et all 120101 ). 10% to the flux density 
(jAguirre et al.ll2010t ). a factor of two uncertainty in the 
opacity (as above), and 20% to the distance. The 68% 
confidence intervals for both dust continuum masses are 
from 80% to 180% of the quoted mass. 

The systematic uncertainty associated with the virial 
mass is difficult to characterize. The basic assumption 
that the clump is a gravitationally bound object in virial 
equilibrium introduces an inherent uncertainty that can- 
not be quantified. Another major uncertainty is the 
virial radius. The true radius of these clumps is in most 
cases unknown, so assigning an uncertainty is somewhat 
questionable. Our typical virial radii are around 0.4 
pc and som e interferometric mil limeter studies of IRDC 
clumps (e.g lBeuther et aLll2005[ ) find sizes of about 0.1 
pc. Therefore, we assign an uncertainty of a factor of 
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4 to the virial radius. The line-width is uncertain by 
about 20% for these data. These uncertainties give us a 
68% confidence interval for the virial masses from 50% 
to 170% of the quoted value. In cases where the adopted 
virial radius is off by a wide margin, the virial mass will 
also be off by that same wide margin (as the calculated 
mass scales linearly with the adopted virial radius). 

The quantifiable systematic uncertainties in the 13 CO 
mass are the excitation temperature, integrated line in- 
tensity, and distance. However, there is a large uncer- 
tainty that we are unable to quantify associated with 
the optical depth, freezeout, excitation conditions, and 
variable abundanc e of 13 CO in cold, dense clumps (e.g. 
iTafalla et~a l. 2002). These will have the effect of lowering 
the 13 CO mass and column density estimates. We assign 
an uncertainty to the excitation temperature of 50% (as 
above), an uncertainty in the integrated intensity of 20%, 
and an uncertainty in the distance of 20%. We find that 
the 68% confidence interval for the 13 CO masses is from 
80% to 120% of the quoted mass, neglecting additional 
uncertainties introduced by optical depth and freezeout 
effects. 

4. ANALYSIS 

4.1. Mass Comparison 

Given the systematic uncertainties, the agreement be- 
tween the various mass tracers is reasonable, as shown 
in Table 2] and Figure O All of the mass tracers increase 
with increasing BGPS 1.1 mm dust continuum mass, ex- 
cept for extinction and 13 CO at the highest BGPS 1.1 
mm masses, as these are associated with UCHII Regions 
(See Table [5J . The virial masses are uniformly higher 
than the BGPS 1.1 mm masses, the ratio of the two 
ranges from 1.1 to 14.4 with a median of 6.4 ± 4.2. A 
detailed comparison of the BGPS 1.1 mm mass with the 
virial mass is given in i j4.3l The MAMBO 1.2 mm masses 
are also uniformly higher than the BGPS 1.1 mm masses, 
though not nearly as high as the virial masses. The ra- 
tio of MAMBO 1.2 mm to BGPS 1.1 mm masses ranges 
from 1.2 to 2.6 with a median of 1.8 ± 0.3. This off- 
set is intrinsic to data taken with MAMBO at 1.2 mm 
vs. Bolocam at 1.1 mm. The reason for this offset is 
being inv estigated by the BG PS team and will be dis- 
cussed in lAguirre etaTl (|2010h . The 13 CO masses agree 
surprisingly well with the BGPS 1.1 mm masses, the ra- 
tio between the two ranges from 0.2 to 1.9 with an me- 
dian of 0.8 ± 0.5. Recall, however, that the 13 CO mass 
is taken over a 46" beam rather than the 33" effective 
beam size of the rest, so we are including more of the 
GMC and the IRDC envelope, and therefore more emis- 
sion. This means that the true comparable 13 CO mass 
would be generally lower than the BGPS 1.1 mm mass. 
The agreement between the two is best at lower BGPS 
1.1 mm masses. At higher BGPS column there is likely 
a higher than assumed BGPS temperature, as well as 
optical de pth effects or freezo ut of 13 CO at high column 
densities ■: I alalia et al.|[2002T) . The extinction mass and 
BGPS 1.1 mm mass agree well, except at the highest 
BGPS 1.1 mm masses (M B gps > 600 M Q ), where the 
sources are associated with UCH II Regions. The ratio 
of extinction to BGPS 1.1 mm mass ranges from 0.1 to 
2.8 with an median of 1.1 ± 0.8. A more detailed com- 
parison of the extinction mass to the BGPS 1.1 mm mass 



is given below. 

4.2. Dust extinction and emission 

The relationship between 8 /im extinction and 1.1 mm 
dust emission is reasonably good, and where it is not, 
there is good indication that warm gas/dust due to a 
hot clump is skewing both extinction and emission mea- 
sures of mass. Recent dust temperature estimates from 
IRathborne et al.l ((2010) corroborate that active clumps 
are indeed warmer than quiescent/intermediate clumps. 
The agreement between 8 fim extinction and 1 . 1 mm dust 
emission masses is poor in the four highest BGPS 1.1 mm 
mass sources, which are associated with UCH II Regions. 
In these regions, the temperature is almost certainly 
warmer than the assumed 15 K, which contributes to the 
very large BGPS fluxes. Additionally, the UCH II Region 
illuminates PAH features in the 8 pum band, boosting the 
8 /im flux and reducing the measured extinction mass. 

Figure [3] shows the extinction column density vs. 
BGPS 1.1 mm column density and we note that while the 
column density derived from extinction and the BGPS 
1.1 mm dust continuum generally agree there are some 
peculiarities. Note that the clumps indicated in Table 
0] are only a small portion of what is plotted in Fig- 
ure |31 Figure [3J plots the entire IRDC, while the table 
only includes the subsections of the cloud chosen for this 
study as clumps. In some sources the extinction flattens 
out as the 1.1 mm column increases, while in others the 
extinction column is higher than that 1.1 mm column. 
These indicate a hot clump (with imperfect temperature 
assumptions for the 1 . 1 mm dust emission and 8 /xm con- 
tamination by PAH illumination) and either a cold clump 
or imperfect foreground/background estimation, respec- 
tively. For the individual clumps, the median of the ratio 
of extinction column to BGPS column is 1.1 ± 0.8 overall, 
0.7 ± 0.9 for active clumps, 1.4 ± 0.6 for intermediate 
clumps and 2.0 ± 0.9 for quiescent clumps. This is in 
agreement with our assessment that high extinction per 
1.1 mm column is characteristic of quiescent clumps, low 
extinction per 1.1 mm column is characteristic of active 
clumps, however, there is still a lot of scatter. Though 
there are some exceptions, in general, dust extinction and 
emission are well correlated for quiescent /intermediate 
clumps, but often diverge for active clumps, particularly 
ones associated with UCHII Regions. 

4.3. A Comparison of Virial and Dust Masses 

Simple calculations show that the measured virial 
masses are uniformly higher than the masses measured 
through other techniques, such as the dust continuum 
emission. The median of the ratio of the virial mass to 
the BGPS 1.1 mm mass is 6.4 ± 4.2, an d ranges from 1 . 1 
to 14.4. Interferometric observations (|Rathbqrne et al.1 
[20071 [2008: W ang et alJl2008tlZhang et al.ll2009L for more 
details on specific sources, please see the Appendix) of 
several of these sources find that most clumps fragment 
into smaller cor es at higher resoluti on, with sizes of about 
0.03 pc or less. iZhang et al.l fj2009f ) find that on the scale 
of ~ 0.02 pc the virial and dust mass estimates agree 
quite well. On the scale of about 1 pc, we find that the 
virial masses are much higher than the dust masses. The 
virial radius required for our measured virial masses to 
match our measured BGPS 1.1 mm masses ranges from 
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Fig. 4. — Ratio of the integrated intensities of N2H+ (J=3-2) to HCO+ (J=3-2). Only pixels greater than 4<r in HCO+ and N2H+ are 
shown. This is a pixel-by-pixel comparison of the HCO" 1 " and N2H + points within the IRDC ellipse. The line drawn has a slope of one 
representing equality. The point in the lower right corner of each plot represents the approximate 1 a error bars for the points in that plot. 



0.02 to 0.25 pc with a mean of 0.09 ± 0.06 pc, simi- 
lar to the interferometric size scales noted above. The 
information presented is consistent with the picture of 
small bound cores (about 0.01 pc in size), surrounded 
by a dense envelope (about 0.1 pc in size) which is also 
bound, engulfed in a diffuse envelope (size scale of about 
1 pc). This dense envelope is where the majority of dense 
gas is emitting and therefore dominates the line width we 
measure, while the diffuse envelope extends out to about 
1 pc and comprises the extended emission we see in the 
dust continuum. In this picture, the line widths we mea- 
sure are primarily from a region that is smaller than the 
virial radius we assign to it, thus explaining the very large 
virial to dust mass ratio. Alternatively, if the material 
producing the large line widths and the dust continuum 
emission are the same, it is possible that something more 
than virialized motion (e.g. magnetic fields, turbulence) 
are playing a role in the line widths or that the clumps 
are not graviationally bound, and possibly ram-pressur e 
confined by converging flows (e.g. iHeitsch et alj l2008f) . 
The question of whether these clouds, clumps, and cores 
are bound is of great importance in understanding cluster 
formation. 



4.4. A Molecular Line Tracer Comparison 

The ratio of HCO + to N2H + in clumps of varying ac- 
tivity levels is of special interest, as it is thought that 
HCO+ will tra ce "hot" clumps whil e N 2 H+ will trace 
"cold" clumps (pOrgensen et alj|200~4f ). HCO+ is created 
by CO in the gas phase, and will also deplete onto dust 
grains at cold temperatures, so it should be more abun- 
dant in hotter environments. N2H + , on the other hand, 
is destroyed by CO in the gas phase, and should be more 
abundant in colder environments. Figure is a pixel- 
by-pixel comparison of N2H+ to HCO + integrated in- 
tensities, and the line is of slope one, representing equal- 
ity. Throughout the majority of the clouds, the ratio of 
HCO + to N2H+ is greater than one. In a few cases (most 
notably G028. 37+00. 07), however, the N 2 H+ integrated 
intensity is greater than the HCO + . This is due to a 
high column of gas and dust, in which HCO + is optically 
thi ck, and self-absorbed which lessens the peak (e.g. Fig- 



15(b)! 



making N2H + brighter by comparison. The 
general trend is that HCO + and N2H + are well correlated 
and that HCO + is slightly brighter. In Figure [5] we see 
that our present sample shows no trend between clump 
activity and the ratio of HCO + to N2H+. It is impor- 
tant to note that we are not resolving individual cores, so 
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Fig. 5. — HCO + vs. N2H + Peak T m t> on a log-log plot. Red is active clumps, green is intermediate and blue is quiescent (See Table[5]l. 
The line drawn has a slope of one, representing equality. Note that there is no noticeable distinction in the ratio of HCO+ to N2H+ for 
active vs. quiescent clumps, but that active clumps are much brighter in HCO + and N2H+ than quiescent clumps. 



adjacent cores with differing chemistry will result in an 
average, rather than resembling distinct core chemistries. 
Additionally, we are comparing line strengths, not abun- 
dances. 

While we see no evidence for a chemical differentiation 
between HCO + and N 2 H + on a ~ 1 pc scale between "ac- 
tive" and "quiescent" clumps, we note that active clumps 
are significantly brighter in HCO + and N 2 H + than inter- 
mediate or quiescent clumps. The increased brightness 
of HCO + and N2H+ in active clumps could be the result 
of an increased column density, higher volume density, 
higher temperature, and/or a larger beam filling factor. 
We argue that an increase in column is not the primary 
effect as not all active clumps have higher dust column, 
yet all active clumps show an increased brightness in 
HCO + and N2H + . The change in filling factor is a pos- 
sible explanation, however the HCO + and N2H" 1 " maps 
show extended emission around all the clumps. Since at 
the typical density of an IRDC (n ~ 10 5 cm -3 ) we are ap- 
proaching the criti cal densities of HCO + and N2H + (see 
iDaniel et al.l 12007b . we suggest that the active clumps 
may have begun to compress and become denser, increas- 
ing the intensity of HCO + and N2H+. This increase 
in density is also responsible for the increased optical 
depth of HCO + in active clumps, as evidenced by the 
self-absorption of HCO + in som e active clumps. Recent 
dust temperature estimates from Rathbor ne et al.l ([2010D 
suggest that active clumps are warmer, which would also 
increase the observed brightness of HCO + and N2H + . 
We observe that HCO + and N2H + are brighter in active 
than quiescent clumps and suggest that this is because 



the active clumps are either warmer, have higher volume 
densities, or both. 

We also compare the HCO + and N2H + peak line in- 
tensities with those of 13 CO in Figure HJ As discussed in 
Section [3~5l while there is good evidence that 13 CO is op- 
tically thick toward many of these clumps, we see that it 
at least traces to some depth in the cold, dense portion 
of the clump. However, the 13 CO is almost certainly 
tracing some part of the GMC and IRDC envelope as 
well. Therefore, we include it as an optically thick com- 
parison point for the dense gas tracers. The 13 CO peak 
varies very little among the different clump types and 
shows no noticeable increase with column, again corrob- 
orating that it is optically thick. The HCO + and N2H+ 
peak temperatures increase somewhat with column den- 
sity, but most notably with activity level. This means 
that the active clumps all have much smaller variation 
in line strength than the quiescent clumps (the line con- 
necting the peak T m f,'s is shorter). The variation in line 
strength among the three tracers can be quantified by the 
difference between the brightest peak line strength of the 
three tracers and the weakest peak line strength of the 
three tracers. The mean of the variation in line strength 
is 4.5 ± 0.8 for active clumps, 4.8 ± 1.1 for intermediate 
clumps, and 5.5 ± 1.8 for quiescent clumps. While the 
numbers are uncertain, the trend is clear, and this is a 
potentially useful tool for identifying active clumps. 

4.5. H II Regions and Stellar Type Limits 

We detect radio continuum sources toward four of the 
clumps. Two of these are unresolved, optically thick 
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TABLE 5 
Clump Activity 



Clump 


3.6 cm Flux a 


H II Region 


"Green Fuzzy" ? b 


24 fimflmt b 


H2O maser? b 


CH3OH maser? b 


Clump c 


Evolutionary* 1 




(mJy) 


Type 


Y/N 


(mJy) 


(Y/N) 


(Y/N) 


Activity 


Stage 


G022. 35+00. 41 


















GLM1 


0.06 


- 


Y 


13 


Y 


Y 


A 


2 


GLM2 


0.36 




N 




N 


N 


Q 


1 


G023. 60+00. 00 


















GLM1(VLA4) 


0.30 


UCH II 


Y 


1058 


Y 


N 


A 


3 


GLM2 


0.12 




N 


198 


Y 


N 


I 


2 


G024. 33+00. 11 


















GLM1(VLA2) 


0.29 


UCH II 


N 


999 


Y 


Y 


A 


3 


GLM2 


0.33 




N 




N 


N 


Q 


1 


G024. 60+00.08 




















n no 

u.uy 




M 
IN 


"1 Q 

lo 


v 
1 


IN 


t 
1 




CJ A/TO 


n no 
u.uy 




v 


..VIZ 


M 
IN 


rN 


T 
1 




G028.23-00.19 


















GLM1 


0.03 




N 




N 


N 


Q 


1 


G028.37+00.07 
















GLM1 


0.18 




N 




N 


N 


Q 


1 


GLM2 


0.12 




Y 


45 


Y 


N 


A 


2 


GLM3 


0.15 




Y 


4 


Y 


N 


A 


2 


GLM4 e 


0.12 




Y 


22 


Y 


Y 


A 


2 


G028.53-00.25 


















GLM1 


0.06 




Y 


16 


N 


N 


I 


2 


G034.43+00.24 


















GLM1(VLA3) 


0.31 


UCH II f 


Y 


1718 


Y 


Y 


A 


3 


GLM2(VLA1) 


10.07 


UCH II 


N 


1401 


Y 


Y 


A 


1 


GLM3 


0.09 




Y 




Y 


Y 


A 


2 



a For point sources this is the peak flux, for the extended source, G034. 43+00. 24: GLM2, the flux is the integrated flux, and for no detection 
(as marked by a — in the 3rd column) it is the 3cr upper limit. 
b IChambers et al.l 02009) 

c A— Active, I— Intermediate, Q— Quiescent 
d Sec J5] for details. 

e G028. 37+00. 07: GLM4 slightly overlaps with the resolved H II region, G028. 37+00. 07 VLA1, see Tablc[F] 
^ G034. 43+00. 24: GLM1 is a marginal detection. It may or may not be an UCH II region. Sec S|4.5l for details. 



TABLE 6 
Radio Continuum Sources 



VLA Source R.A. DECL. Peak 3.6 cm Flux Total 3.6 cm Flux cr Point Source a 

(J2000) (J2000) (mjy/bcam) (mJy) (mjy/beam) (Y/N) 

G022.35+00.41 



G023.60+00.00 



VLA1 


18:34:33.1 


-08:15:26.8 


24.3 


28.4 


0.3 


Y 


VLA2 


18:34:09.6 


-08:17:49.4 


1.5 


68.9 


0.1 


N 


VLA3 


18:34:12.4 


-08:19.01.5 


1.1 


1.3 


0.1 


Y 


VLA4 (GLM1) 


18:34:21.1 


-08:18:12.4 


0.3 


0.3 


0.1 


Y 


G024.33+00.11 














VLA1 


18:35:24.0 


-07:37:37.9 


2.8 


3.6 


0.5 


Y 


VLA2 (GLM1) 


18:35:08.1 


-07:35:04.0 


0.29 


0.4 


0.02 


Y 


G024.60+00.08 














G028.23-00.19 














VLA1 


18:43:33.7 


-04:11:48.2 


0.24 


0.4 


0.03 


Y 


G028.37+00.07 














VLA1 


18:42:52.9 


-04:00:09.2 


3.0 


49.9 


0.5 


N 


VLA2 


18:42:41.4 


-04:02:16.9 


1.7 


9.6 


0.4 


N 


VLA3 


18:43:03.5 


-04:00:13.2 


1.98 


2.2 


0.08 


Y 


VLA4 


18:42:37.1 


-04:02:02.3 


1.5 


1.7 


0.2 


Y 


VLA5 


18:43:04.5 


-04:00:15.0 


1.16 


1.5 


0.08 


Y 


VLA6 


18:42:52.4 


-03:59:07.0 


0.38 


0.1 


0.04 


Y 


G028.53-00.25 














G034.43+00.24 














VLA1 (GLM2) 


18:53:18.7 


+01:24:47.2 


8.79 


10.1 


0.07 


N 


VLA2 


18:53:08.4 


+01:29:33.9 


1.00 


1.1 


0.08 


Y 


VLA3 (GLM1) 


18:53:18.1 


+01:25:25.4 


0.31 


0.4 


0.06 


Y 



a Extended sources are plotted in Figurcs |13|26l 
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Fig. 6. — A comparison of the peak main beam temperatures for 13 CO (cross), HCO + (triangle), and N2H+ (diamond) vs. BGPS column 
density. Each dashed line represents a single clump, and the color of the line denotes whether it is active (red), intermediate (green), or 
quiescent (blue). Note that 13 CO is relatively flat with increasing column density, corroborating that it is optically thick. The variation 
in line strength between the various tracers (the length of the dashed line drawn) is much smaller for active vs. quiescent clumps and is a 
potentially useful diagnostic for understanding star formation activity. 
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UCH II regions. One is a slightly resolved, optically 
thin UCH II region. The final source is an unresolved, 
marginal detection which may be an UCH II region. The 
peak fluxes of these regions are given in TableEl as well as 
the limiting (3cr) flux in the clumps without detections. 
The properties of all the detections (outside or inside the 
clumps) are reported in Table [5] The sources are named 
based on the IRDC in which they were found (though 
they are not necessarily associated with the cloud) and 
in order of decreasing brightness (VLAl, VLA2, etc.). 
All the Gaussian-type sources were fit with 2-D ellipti- 
cal Gaussians, and the fluxes reported are from the fits. 
The a reported is the r.m.s. of the fit residuals. For the 
sources with non-Gaussian morphology, the fluxes were 
measured using a circular aperture and an outer annulus 
of equal area. The mean of the flux in the outer annulus 
is assumed to be the background level and is subtracted 
from the source flux, while the standard deviation in the 
outer annulus is our quoted a. Some of these detec- 
tions are more evolved H II regions, peripheral to the 
IRDC. Others are point sources unassociated with any 
GLIMPSE emission. These sources may be unresolved 
H II regions that happen to be in our field of view but 
unassociated with the mid-IR, or they may be other types 
of Galactic or extragalactic radio continuum sources. We 
restrict our analysis to the four radio continuum sources 
associated with the millimeter clump peaks, as these are 
almost certainly UCH II regions in the IRDC. 

The UCH II region radio continuum flux is directly 
related to the number of Lyman continuum photons 
(Q) of the ionizing star if the region is optically thin. 
Through comparison with the 20 cm and 6 cm ra- 
dio continuum data from the Multi- Array Gala ctic 
Plane Imaging Survey (MAGPIS; iWhite et al.ll2005h we 
find that G034.43+00.24: GLM2 (VLAl) is slightly re- 
solved (D ~ 0.14 pc) and optically thin at 3.6 cm. 
G034.43+00.24: GLM1 (VLA3) is a marg inal detection in 
both o ur survey and the 6 cm survey by FShcphcrd et al.l 
(2004). The measured flux densities at 3.6 and 6 cm 
are not consistent with the source being an UCH II re- 
gion. However, the fluxes reported are highly uncer- 
tain, and systematic errors could be large enough to ac- 
count for this inconsistency if it is in fact an UCH II 
region. Since an optically thick UCH II region is the 
most likely explanation for a radio continuum source as- 
sociated with star formation in an IRDC clump, we as- 
sume that this is the explanation for the radio contin- 
uum source. We use the 3.6 cm flux to calculate the 
Lyman continuum flux of an unresolved, optically thick, 
UCH II region, but acknowledge that this assumption 
may be incorrect. For G024.33+00.11: GLM1 (VLA2) 
and G023.60+00.00: GLM1 (VLA4) the 6 cm upper lim- 
its are consistent with the source being an optically thick 
UCH II region at 3.6 cm. 

In the optically thin case (G034.43+00.24: GLM2 
(VLAl)), we fit a thermal bremsstrahlung curve to the 
3.6, 6, and 20 cm points. We find an emission mea- 
sure EM = 4.4 x 10 6 cm -6 pc, a turnover frequency 
v{j = 1) = 1.3 GHz, and Q = 6.1 x 10 45 Lyman contin- 
uum photons per s econd. This translate s to a B0. 5 star 
( Vacca et al.l Il996^. in agreement with Shepherd et al.l 
20041) . 

For the unresolved, optically thick sources, we first de- 



rive a source radius 



S„c 2 
2v 2 kT P 



-4D Z 



1/2 



(20) 



where S„ is the radio continuum flux, v is the band-center 
frequency, T e is the ionized gas temperature, and D is 
the distance. We assume a T e = 8000 K ionized gas, 
typical of an H II region. We conservatively assume that 
the observed frequency is the turnoff frequency (where r 
= 1). This gives us a lower limit for the emission mea- 
sure (EM), electron density (n e ), and number of Lyman 
continuum photons (Q). If our observing frequency is the 
turnover frequency (v) then r = 1 in our expres sion for 
the emission measure (IWood fc ChurchwelHll989l) 

EM(cm~ 6 pc) 



8.235 x W- 2 a(v,T e )Te 1S5 v- 21 ' 

( 21 ) 

The lower limit for the number density of electrons, 
which we take to be equal the number density of ions, 
is 



EM 
~2r r ' 



(22) 



The number of Lyman continuum photons is given by 
the Stromgren sphere equation, 



Q = -irr a B n e 



(23) 



where as is the recombination rate coefficient, 3.1 
x 10~ 13 cm 3 s" 1 for T e = 8000 K. The Stromgren 
sphere equation assumes that the H II region is spher- 
ical, in equilibrium, that there is one electron per 
ion, and that each energetic photon ionizes an atom 
or molecule. Since the observed structure is clumpy, 
the Stromgren sphere approximation is probably not a 
very good one. Using this method for the two (pos- 
sibly three) optically thick, unresolved UCH II regions 
yields a l og(Q) of 46 . 14, 4 6.55, and 46.15 and stel- 
lar types (jVacca et al.l fl99l of B0, 09.5, and B0 for 
G023.60+00.00: GLM1 (VLA4), G024.33+00.11: GLM1 
(VLA2), and G034.43+00.24: GLM1 (VLA3) respec- 
tively. 

4.6. Clump Activity 

A major driving question in the study of IRDCs is the 
evolutionary sequence of the clumps and clouds: identi- 
fying the different evolutionary stages and their relative 
lifetimes. This requires an understanding of the star for- 
mation activity. We emply four different star formation 
tracers: 1) an embedded 24 /jm point source, 2) "green 
fuzzies," regions of enhanced, extended 4.5 /jm emission, 
indicative of shocks and outflows, 3) H2O and Class I 
CH3OH maser emission, and finally, 4) Ultra-Compact 
(UC) H II regions, an unambiguous indication that a 
massive star that has "turned on." Tracers of star for- 
mation activity in these clumps are summarized in Table 

El 

Each of these tracers is sensitive to somewhat differ- 
ent environments. A 24 fim point source traces material 
accreting onto forming stars, the gravitational contrac- 
tion of a young stellar object (YSO), or even the nuclear 
energy emitted once the star has "turned on." "Green 
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Fig. 7. — 24 fj,m flux (as measured by Chambers ct al. 2009) vs. dust temperature (from Rathbomc ct al. 2010) for the red: "red" 
clumps, orange: clumps with three star formation tracers, green: two star formation tracers, and blue: one star formation tracer. Note the 
that "red" clumps are associated with brighter 24 (im flux and that in general, the more tracers of star formation activity, the warmer the 
clump. 



fuzzies" (jChambers et all 120091: ICveanowski et al.l 120081 
also called Extended Green Objects, EGOs) are thought 
to arise from shock-excited spectral line e mission 
(|Marston et al.ll2004j: iNoriega-Crespo et al.112001 due to 
outflows from young stars. Man y of the HCO + and 
N2H + spectra (see Figures 11(b) to 26(d) I show asym- 



metric line profiles commonly associated with outflows, 
many of which (e.g. G034.43+00.24 GLM3, Figure 



26(d) I are associated with "green fuzzies," further impli- 
cating "green fuzzies" as positive outflow tracers. H2O 
(22.23 GHz) and CH 3 OH (Class I 24.96 GHz) masers 
are well-known signposts of star formation activity, likely 
formed in shocks and outflows. Finally, we include the 
presence of 3.6 cm radio continuum emission which when 
associated with a millimeter peak or mid-IR emission is 
a clear indicator of an UCH II region. 

Together, these tracers are sensitive to warm, embed- 
ded dust indicative of accretion, shocks and outflows, 
and thermal bremsstrahlung from newly formed H II re- 



gions. While these are all excellent tracers of star for- 
mation, none are perfect, so we would not expect to see 
every tracer in every actively star-forming region. View- 
ing angle or dust obscuration could play a significant 
role in a non-detection. We assign the designation "ac- 
tive" to clumps that exhibit three or four signs of ac- 
tive star formation ("green fuzzy," 24 /xm point source, 
UCH II region, or maser emission). This ensures that 
each "active" clump has either an H II region (unam- 
biguously star-forming) or at least two outflow tracers 
(maser emission or "green fuzzy") and a 24 /im point 
source. We assign the designation "quiescent" to clumps 
that exhibit no signs of active star formation. We re- 
serve the title "intermediate" for clumps which exhibit 
one or two signs of active star formation. Intermediate 
clumps, therefore, may have only shock/outflow signa- 
tures or a 24 /im point source. These categorizations 
represent a crude separation between the quiet, massive- 
starless clumps and the active, star-forming clumps, not 
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an evolutionary sequence. Evolutionary sequences are 
discussed in <J5] 

The detection of several UCH II regions embed- 
ded in IRDC clumps unambiguously shows that some 
IRDCs are forming massive stars. We note that all the 
clumps associated with UCH II regions have significantly 
brighter (> 1 Jy) 24 fim flux. All but one of the UCH II 
regions are associated with a CH3OH maser and two of 
the four have a "green fuzzy." The lack of "green fuzzies" 
and a CH3OH maser in some sources may indicate that 
the outflow stage has ceased, or could simply be related 
to an unfavorable viewing angle. We also note that the 
brighest UCH II regions (G034.43+00.24: GLM2 and 
near the edge of G028.37+00.07: GLM4) are associated 
with "diffuse red clumps," regions of extended, enhanced 
8 /xm emission. The ionizing radiation from a young H II 
region will excite PAH features in the 8 /im band. As the 
H II region grows it will begin to clear some of the sur- 
rounding dust and gas and reveal itself as a "diffuse red 
clump." As the H II region continues to evolve, there will 
be PAH destruction by the strong UV flux and it will no 
longer have enhanced 8 /xm emission. The clumps which 
contain UCH II regions, but are not classified as "red" 
are either in an earlier evolutionary phase or are perhaps 
still obscured by the dust at 8 /xm. There are no "dif- 
fuse red clumps" in our sample which are not associated 
with UCH II regions, so these are potentially very useful 
tr acers of young H I I regi ons. 

IChambers et al.1 (|2009l ) discussed the use of "red 
clumps" (regions of enhanced 8 /xm emission) as trac- 
ers of embedded H II regions. Since an enhancement 
at 8 /im could presumably be caused by any significant 
UV flux and most high column density clumps are nearly 
optically thick at 8 /im, we suggest that only the more 
diffuse "red clumps" are necessarily associated with H II 
regions. Only when the forming star has cleared out some 
dust obscuration and is no longer associated with a mil- 
limeter peak can it be identified as a "diffuse red clump." 
Other "red clumps" could be embedded UCH II regions 
that have not cleared out their envelopes, or they could 



be PAH nebulae excited by main sequence or evolved AB 
stars. 

Figure [7] shows th e 24 /im flux (as measured by 
IChambers et al.l [2009D vs. the dust temperature (mea- 
sured through SEP fit s in the sample of 100 clumps by 
iRathborne et al.l I201(jh . Notice that the "red" clumps 
are associated with brighter 24 /xm flux and that in gen- 
eral, the more star formation tracers, the warmer the 
clump is. This general trend indicates that these tracers 
combined are successfully measuring star formation ac- 
tivity. However, the large scatter may be due to viewing 
angle biases, mis-assigned star formation tracers, or the 
inadequacy of single-temperature SED fits. 

5. DISSCUSSION OF AN EVOLUTIONARY SEQUENCE FOR 
IRDC CLUMPS 

Based on a samp le of 190 IRDC clumps, 
IChambers et alj (|2009f ) proposed an evolutionary 
sequence for IRDC clumps which begins with a quies- 
cent clump (in their definition, this means no "green 
fuzzy" or 24 /im point source), transitions into an 
active clump (the clump contains a "green fuzzy" and 
a 24 /xm point source), and finally, becomes a "red" 
(enhanced 8 /xm emission) clump. Our radio continuum 
data support this sequence, especially that the final 
evolutionary stage is a "red" clump, in which the 
H II region is exciting the PAH spectral feature in 
the 8 /xm band. However, we suggest that the term 
"red" clump be reserved for more diffuse, evolved "red" 
clumps, not associated with a millimeter peak. This 
prevents any confusion with embedded B or A stars 
which have UV luminosity signficant enough to excite 
the PAH emission in the 8 /xm band, but not enough 
to produce an UCH II region. Some clumps with H II 
regions are "red" and some are not, which may indicate 
different environments or evolutionary stages. Clumps 
with UCH II regions tend to have a significantly brighter 
24 /xm source flux (> 1 Jy). 

Presented in Figure [8] is a cartoon of star formation 
tracers and how they may possibly relate to evolution- 
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Fig. 9. — Left: For each star formation activity tracer used, we plot a dust temperature (as measured by Rathborne ct al. 2010) histogram 
of the clumps with ( solid line) and without (dashed line) that star formation tracer (as measured byJChambcrs ct al. 2009) from the sample 
of lRath borne et al. (201(J). Note that all the star formation tracers exist over a wide range of temperatures, except for "red" clumps which 
are found primarily at higher temperatures. The top four panels are Stage 2 sources while the "red" clump panel are Stage 4 sour c es (se e 
Section[5]l. Right: A dust temperature histogram of possible evolutionary stages in the sample of 100 clumps from [Rathborne ct al. (2010). 
This sample is likely incomplete at the low-temperature end, as many of the coldest sources would also be faint and below the detection 
threshhold for the SED fits. 



ary stages. Stage 1 is a quiescent clump; a cold, millime- 
ter peak with no signs of active star formation. Stage 
1 can be split into two sub-categories: a) a quiescent 
clump which has not yet begun to collapse and b) a 
quiescent clump that shows collapse signatures, but no 
signs of active star formation. The next stage is a hot- 
ter, denser clump that shows signs of active star forma- 
tion. Stage 2 is a clump with at least 1 sign of active 
star formation, a 24 /im point source, "green fuzzy," or 
a H2O or CH3OH maser. We cannot infer any sequen- 
tial nature between these observational signatures at this 
time. However, the authors hypothesize that, in analogy 
with low- mass star forming regions, an outflow ("green 
fuzzy" ) and 24 fim point source will occur first (the one 
seen first depends on viewing angle), followed by an H2O 
maser, and then a CH3OH maser as the accreting proto- 
stars gain mass. 

The third stage in this sequence is the ignition of an 
UCH II region (see Figure [5]). A Stage 3 source will 
have a strong (> 1 Jy) 24 /im point source and possibly 
mas er emission and outflow indicators. Modern theory 
(e.g. lMcKee fc Ostr ikcr 2007) suggests that massive stars 
reach the main sequence while still accreting, because the 
accretion time scale is greater than the Kelvin- Helmholtz 
timescale. Thus, a main-sequence B star may continue 
to accrete on the main sequence, eventually becoming an 
O star. A Stage 3 source, therefore, may remain compact 
and continue accreting for some time. The fourth stage 
in this evolutionary sequence is the expansion of the H II 
region, producing a diffuse "red" clump (see Figure H]). 
However, this is not a necessary condition for massive 
star formation, as in the densest cores, the H II region 
may be quenched by the remaining envelope (|Ketoll2003f) 
for some time and a "diffuse red clump" will not form. 
For UCH II regions not trapped by their surrounding 
envelope, the UCH II region will expand as the radiation 
field increases or instabilities develop. This expansion 
evacuates a dust cavity, creating a diffuse enhancement 
at 8 /im, that is spatially distinct from any millimeter 



continuum peaks. In Table [5] we assign each clump to 
these coarse evolutionary stages. 

One potential test of the sequential nature of star for- 
mation tracers is the clump temperature. We suggest 
that as a clump evolves, its temperature monotonically 
increases with time. In a statistical sense, warmer ob- 
jects should be more evolved than col der objects. We 
use th e dust temperatures derived by IRathborne et al.l 
(2010) in a sample of 10 clumps and the star f ormation 
tracers as catalogued bv lChambers et ail (|2009fl to inves- 
tigate how dust temperatures change with star formation 
tracers in Figure |H1 All of the individual star formation 
tracers exist over a surprisingly wide range of tempera- 
tures, except for "red" clumps which are found primarily 
at higher temperatures. Figure [7] shows that the more 
star formation tracers a clump has, the warmer it is on 
average, however, no sequential nature can be inferred 
between the individual star formation tracers in Stage 2. 
The lower panel of Figure [3] shows that the broad catego- 
rizations between Stage 1 (quiescent), Stage 2 (interme- 
diate/active), and Stage 4 ("diffuse red clump") clumps 
are well-separated by temperature, and likely represent 
a true evolutionary sequence (we note that many Stage 1 
clumps were too faint to reliably fit an SED, so the true 
histogram is likely even further to the colder side of the 
histogram) . These data support the coarse separation of 
Stage 1, Stage 2, and Stage 4 clumps as an evolutionary 
sequence, however, no further sequential nature of star 
formation tracers can be inferred. 

When discussing evolutionary sequences, there is al- 
ways the added uncertainty of whether we are truly see- 
ing an evolutionary sequence or just different environ- 
ments. We point out that it is likely that not all Stage 
1 (quiescent) clumps will form Stage 2 (active) clumps, 
however, we are suggesting that all Stage 2 clumps form 
from Stage 1 clumps and so on. Some quiescent clumps 
may be forming low/intermediate mass stars which will 
remain undetected by our star formation tracers. Alter- 
natively, some quiescent clumps may be transient density 
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enhancements which will eventually disperse. This makes 
the question of lifetimes more complicated. The lack of 
a sequence between the observed star formation tracers 
is a bit surprising, however, the fact that more star for- 
mation tracers indicates a higher temperature indicates 
that viewing angle/obscuration may be playing a large 
role in the observed scatter. Perhaps only when a clump 
is more evolved and somewhat less obscured by dust can 
one detect many of these star formation tracers. Higher 
sensitivity and better spatial resolution observations, as 
well as a better control of systematics may reveal some 
sequential nature between the observed star formation 
tracers. We note that th e ove rlap of the Stage 2 and 
Stage 4 sources in Figure 9(b) hint at some number of 
UCH II regions undected as "red" clumps or some num- 
ber of "red" clumps incorrectly associated with UCH II 
regions. 

6. CONCLUSION 

We examined a sample of 17 clumps within 8 IRDCs 
using existing infrared, millimeter, and radio data, as 
well as new VLA radio continuum and HHT dense gas 
spectral data. Our main conclusions are: 

• 8 /im extinction mass and BGPS 1.1 mm mass are 
complementary tracers of mass in IRDCs, except 
for the most active clumps (notably those contain- 
ing UCH II regions), for which both mass tracers 
suffer biases. 

• The measured virial masses of IRDC clumps are 
uniformly larger than the dust continuum masses 
on the scale of ~ 1 pc. 

• We do not detect a chemical differentiation of 
HCO + and N2H + between active and quiescent 
clumps on the scale of ~ 1 pc. However, both 
HCO + and N2H + are brighter in active clumps, 
due to an increase in temperature and/or density. 

• We report the identification of four UCH II re- 
gions embedded within IRDC clumps. We find that 
UCH II regions are associated with bright (> 1 Jy) 



24 /an point sources, and the brightest UCH II re- 
gions are associated with "diffuse red clumps" (an 
extended enhancement at 8 /im). 

• We discuss an evolutionary sequence for cluster- 
forming clumps and find that the sequential sepa- 
ration of the broad evolutionary stages, quiescent, 
active, and "diffuse red clumps," is supported by 
dust temperatures, however, no sequential nature 
can be inferred between the individual star forma- 
tion tracers. 
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APPENDIX 
INDIVIDUAL CLUMPS AND CLOUDS 

This appendix explores the properties of each cloud in more detail. Presented are GLIMPSE three-color images of 
each cloud with MAMBO 1.2 mm contours and MIPS GAL 24 ^im images with BGPS 1.1 mm contours. For clouds 
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Fig. 10.— G022. 35+00.41: Left: GLIMPSE three-color image, red is 8 /mi, green is 4.5 /im and blue is 3.6 /jm with MAMBO 1.2 mm 
contours overlaid. The contours are on a log scale from 15 to 410 mjy beam' 1 . Right: MIPSGAL 24 fim image with BGPS 1.1 mm 
contours. The contours are on a log scale from 0.12 to 0.7 Jy beam 
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Fig. 11. — G022. 35+00. 41: Left: GLIMPSE 8 fim overplotted with the BGPS beam-sized apertures that were used to determine clump 
masses. Right: HCO+, N 2 H+ and 13 CO spectra in clumps GLM1 (Stage 2) and GLM2 (Stage 1). 




Fig. 12.— G023. 60+00. 00: Left: GLIMPSE three-color image, red is 8 /im, green is 4.5 /im and blue is 3.6 fj,m with MAMBO 1.2 mm 
contours overlaid. The contours are on a log scale from 30 to 440 mjy beam- 1 . Right: MIPSGAL 24 fim image with BGPS 1.1 mm 
contours. The contours are on a log scale from 0.12 to 1.1 Jy beam -1 . 
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Fig. 13. — G023. 60+00. 00: Top Left: BGPS 1.1 mm continuum dust emission overlaid with linear HCO + contours from 1 to 10 K km s _1 . 
Top right: BGPS 1.1 mm continuum dust emission overlaid with linear N2H+ contours from 1.2 to 8.2 K km s — 1 Middle: GLIMPSE 8 
overplotted with red BGPS beam-sized apertures that were used to determine clump masses. The blue X's are VLA 3.6 cm point sources, 
and extended 3.6 cm sources are depicted as blue contours on the adjacent box. Right: HCO + , N2H + and 13 CO spectra in clumps GLM1 
(Stage 3) and GLM2 (Stage 2). 
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Fig. 14.— G024. 33+00. 11: Left: GLIMPSE three-color image, red is 8 /im, green is 4.5 /im and blue is 3.6 fj,m with MAMBO 1.2 mm 
contours overlaid. The contours are on a log scale from 30 to 1280 mjy beam -1 . Right: MIPSGAL 24 fim image with BGPS 1.1 mm 
contours. The contours are on a log scale from 0.12 to 2.5 Jy beam -1 . 

where HCO + and N2H+ maps were made, they are presented as contours on a BGPS 1.1 mm greyscale image of the 
cloud. Also presented are the position and sizes of the apertures used for clump masses, the location of VLA 3.6 cm 
point sources, and where resolved, contours of 3.6 cm emission on top of a GLIMPSE 8 /im image. We also present 
13 CO, HCO + , and N2H + spectra at each clump position. 

G022.35+ 0041 

G022. 35+00. 41 (see Figures [T0l and [TTj) is actually composed of two clumps along the line-of-sight, one at 3.6 kpc 
and the other at 4.8 kpc. At least one of the two clumps is truly an active star-forming region. G022. 35+00. 41: GLM1 
shows a bright millimeter peak, masers, a "green fuzzy," an d a 24 /xm point source. Interferometric measurements at 1 
and 3 mm on the IRAM Plateau de Bure Interferometer by Rathbo rne et al.l (|2007f) show that GLM1 contains 2 cores, 
less than 0.026 pc each. The second clump (and associated filament) is quiescent and has diffuse millimeter emission. 
While the superposition of two clouds along the line-of-sight lessens the intrinsic column of either, this example shows 
us that in the confused inner Galaxy, it is not terribly uncommon to have two spatially coincident dense chimps. 



G023. 60+00.00 

G023. 60+00. 00 (see Figures [12] and fT3]) is a particularly interesting example with regard to the comparison of 
8 /im extinction and dust emission. The extinction masses of the two clumps are very close (100 and 120 M ), 
however, the BGPS 1.1 mm masses differ by almost a factor of two (140 and 80 M©). GLM1 is an active clump, 
so it is likely warmer and denser than GLM2, an intermediate clump. Also, there is a bright H II region to the 
southwest corner of the image shown in Figure [12] which increases the foregound 8 /im emission. Al so, the HC O + an d 
N2H + maps show only diffuse emission around GLM2 and a bright peak toward GLM1 (see Figures [l3(a)1 and p(b)t , 
tracing the dense, hot gas associated with an active clump. 

At 2.9 x 1.8 pc (major x minor axes of the IRDC ellipse, see Figured]), G023. 60+00. 00 is on the smaller side 
for an IRDC, and shows a filamentary morphology. It is near two evolved H II regions an d a co mplicated web of 
mid-infrared bubbles, filamentary IRDCs, and Photo-Dissociation Regions (PDRs; see Fig ure 13(c) ). The loca tion of 
G023. 60+00. 00 is at 3.6 kpc, right amidst the active star formation in the molecular ring (jjackson et al. 2006). 
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Fig. 15. — G024. 33+00. 11: Left: GLIMPSE 8 fim overplotted with red BGPS beam-sized apertures that were used to determine clump 
masses. The blue X's are VLA 3.6 cm point sources. Right: HCO+, N 2 H+ and 13 CO spectra in clumps GLMl (Stage 3) and GLM2 (Stage 
1). The HCO + and N2H+ spectra have been multiplied by a factor of two, in order to see them more clearly relative to the 13 CO. 




Fig. 16.— G024. 60+00. 08: Left: GLIMPSE three-color image, red is 8 /im, green is 4.5 /im and blue is 3.6 /an with MAMBO 1.2 mm 
contours overlaid. The contours are on a log scale from 15 to 360 mjy beam -1 . Right: MIPSGAL 24 fira image with BGPS 1.1 mm 
contours. The contours are on a log scale from 0.12 to 0.8 Jy beam -1 . 
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Fig. 17. — G024. 60+00. 08: two distinct velocity components/clumps along the line-of-sight: Top left: BGPS 1.1 mm continuum dust 
emission overlaid with linear HCO + contours from 1.2 to 4.4 K km s from GLM1. Top right: BGPS 1.1 mm continuum dust emission 
overlaid with linear N2H+ contours from 1.5 to 5 K km s — 1 from GLM1. Bottom left: BGPS 1.1 mm continuum dust emission overlaid 
with linear HCO + contours from 1.8 to 9.8 K km s from GLM2. Bottom right: BGPS 1.1 mm continuum dust emission overlaid with 
linear N2H + contours from 1.5 to 5.5 K km s — ^rom GLM2 
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Fig. 18. — G024. 60+00. 08: Left: GLIMPSE 8 fira overplotted with red BGPS beam-sized apertures that were used to determine clump 
masses. The blue X's are VLA 3.6 cm point sources, and extended 3.6 cm sources are depicted as blue contours in the adjacent box. Right: 
HCO+, N 2 H+ and 13 CO spectra in clumps GLM1 (Stage 2) and GLM2 (Stage 2). The HCO+ and N 2 H+ spectra have been multiplied 
by a factor of two, in order to see them more clearly relative to the 13 CO. Note the distinctly different clump velocities 
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Fig. 19.— G028.23-00.19: Left: GLIMPSE three-color image, red is 8 /im, green is 4.5 /im and blue is 3.6 jum with MAMBO 1.2 mm 
contours overlaid. The contours are on a log scale from 30 to 112 mjy beam -1 . Right: MIPSGAL 24 fim image with BGPS 1.1 mm 
contours. The contours are on a log scale from 0.12 to 0.64 Jy beam -1 . 



G024. 33+00. 11 

G024. 33+00. 11 (see Figures IT4l and fTS"]) contains a dense, active clump and more diffuse filamentary structure. There 
is a ridge between dark 8 /im extinction and bright 8 /im emission at the southern edge of the cloud. At a distance of 
5.9 kpc, this cloud is among the most distant in our sample. G024. 33+00. 11 was not mapped in HCO + and N2H+, 
but it would be a good candidate for a future mapping project, to see if the dense ridge is associated with the IRDC. 
Additionally, portions of the cloud are dark at 8 /im, but show very little millimeter emission. 

G024. 33+00. 11: GLM1 is a particularly good example of a dense, active clump. Interferometric measurements at 1 
and 3 mm on the IRAM Plateau de Bure Inter ferome ter bv lRathborne et all (|2007t ) show that G024. 33+00. 11: GLM1 
has a single core smaller than 0.035 pc. Figure 15(b) exemplifies the increase of HCO + and N2H 4 " in an active clump, 
and the self-absorption of 13 CO and HCO + in a particularly dense environment. This is in contrast to the weak 
HCO+ and N 2 H+ detections in G024.33+00.11: GLM2. 



G024. 60+00. 08 

G024. 60+00. 08 (see Figures (THl El and lTSl) is the second (of two) examples in our sample which is comprised of two 
unique clumps along the line-of-sight, rather than one contiguous object. The two clumps are spatially adjacent with 
only a slight overlap of diffuse emission, which is an astonishing coincidence. The HCO + and N2H + maps confirm (see 
Figure fl7|) this chance alignment. 

At 6.0 kpc, G024. 60+00. 08: GLM1 is the most distant clump in our sample. The 8 /jm extinction does not appear 
abnormally large, but given the amount of foreground emission between us and the clump, it has one of the highest 
extinction masses in our sample. G024. 60+00. 08: GLM2, on the other hand, is located in the molecular ring at 3.4 
kpc. This clump shows self-absorption of HCO + , and an asymmetric line profile in 13 CO, indicative of outflows. These 
clumps are both catego rized as intermed i ate. I nterferometric measurements at 1 and 3 mm on the IRAM Plateau de 
Bure Interferometer bv lRathborne et all (|2007t ) show that G024. 60+00. 08: GLM1 contains 5 cores, less than 0.041 pc 
each, while GLM2 contains 3 cores, smaller than 0.024 pc each. 



G 028. 23-00. 19 

G028. 23-00. 19: GLM1 (see Figures fl9l and |20|) is the prime example of a starless IRDC clump. It is typical of an 
IRDC clump in size, mass, 8 /im extinction, and in having a relatively compact millimeter core. However, this clump 
has absolutely no signs of star formation. The detections of HCO + and N2H + in this source are incredibly weak. The 
dust column is comparable to many other clumps with bright HCO + and N2H + , the primary difference being the lack 
of density and a heating source in G028. 23-00. 19: GLM1. The 8 /im point source just north of GLM1 may be an 
indication of nearby, recent star formation. Given the compact nature of GLM1, its high column density, and nearby 
star formation, it seems likely that it is a quiescent precursor to massive star formation. 



G028. 37+00. 07 

G028. 37+00. 07 (see Figures |2"T1 and |2"2"|) is an extremely varied cloud. It was identified bv lSimon et al.l (|2006a|) as the 
darkest (by contrast) cloud in the First Galactic quadrant. Due to a particularly bright background around the cloud 
at 8 /im (as discussed in ^3.2j) the extinction column density is over-estimated as compared with the BGPS 1.1 mm 
column density. This cloud is the most massive (by all tracers) in our sample. G028. 37+00. 07 is composed of dense 
filaments and surrounded by an extremely bright mid-IR background, including nearb y evolved H II regions. 



G028. 37+00. 07 hosts a prime example of a "diffuse red clump." In Figure [21(a)] we see a "diffuse red clump,'' 
spatially separated from the bright millimter peak. The H II region producing this red" clump is asymmetric. Away 
from the millimeter clump, the H II region has expanded and become more diffuse and toward the millimeter clump 
the H II region remains confined by the high-density. 
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Fig. 20. — G028. 23-00. 19: Left: GLIMPSE 8 /tm overplotted with red BGPS beam-sized aperture that was used to determine the clump 
mass. The blue X's are VLA 3.6 cm point sources. Right: HCO+, N 2 H+ and 13 CO spectra in clump GLM1 (Stage 1). The HCO+ and 
N2H + spectra have been multiplied by a factor of ten, in order to see them more clearly relative to the 13 CO. 
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Fig. 21.— G028. 37+00.07: Left: GLIMPSE three-color image, red is 8 /mi, green is 4.5 /im and blue is 3.6 /an with MAMBO 1.2 mm 
contours overlaid. The contours are on a log scale from 30 to 1340 mjy beam -1 . Right: MIPSGAL 24 fim image with BGPS 1.1 mm 
contours. The contours are on a log scale from 0.2 to 2.2 Jy beam -1 . 

The structure of HCO+ and N 2 H+ in G028.37+00.07 is a great example of the differing ratio of HCO+/N 2 H+. The 
HCO + and N2H 4 " maps show a similar overall structure, but the intensities of the two tracers vary across different 
clumps. In HCO + , GLM2 is the brightest, and in N2H+, GLM4 is by far the brightest. As pointed out previously, 
however, this is not due to chemical differentiation. Rather, the HCO + and even the N2H 4 " become optically thick 
(and self-absorbed) in the densest clump, GLM4. The self-absorption of HCO + makes N 2 H + appear brighter by 
c omparison. 

iWang et all ([20011 and lZhang et~aH (|2009l ) observed G028.37+00.07: GLM2 and GLM4 (referred to as PI and P2) 
with the VLA and SMA, respectively. In GLM2, they find 5 cores along the filament with line widths of about 1.2 
km s -1 and temperature of about 15 K. In GLM4, they find 2 cores with line widths of about 4.3 km s _1 , temperatures 
of about 30 K, and a rich molecular spectra. 

G028.53-00.25 

G028. 53-00. 25 (see Figures [23] and [24]) is a curious example of a source with a very small extinction column as 
compared with the BGPS 1.1 mm column. This cloud either has an exceptionally low background or high foreground. 
The region surrounding G028. 53-00. 25 is mostly devoid of extended H II regions, PDRs, and other strong emitters of 
diffuse 8 /zm emission. Due to the paucity of nearby diffuse 8 emission, the background estimate is decreased, and 
the extinction column estimate too low. 

G028. 53-00. 25 is a fairly quiescent cloud with few indicators of star formation. The clump, GLM1, is identified as 
intermediate due to the presence of a "green fuzzy" (though a fairly compact, faint one) and 24 /im point source. GLM1 
shows very weak HCO + and N2H 4 " emission (the OTF-maps suffer from low signal to noise). G028. 53-00. 25: GLM1 
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Fig. 22. — G028. 37+00. 07: Top left: BGPS 1.1 mm continuum dust emission overlaid with linear HCO+ contours from 1.2 to 9 K km s _1 . 
Top right: BGPS 1.1 mm continuum dust emission overlaid with linear N2H+ contours from 0.8 to 10 K km s — 1 . Middle: GLIMPSE 8 /im 
overplotted with the BGPS beam-sized apertures that were used to determine clump masses. The blue X's are VLA 3.6 cm point sources, 
and extended 3.6 cm sources are depicted as blue contours in the adjacent boxes. Right: HCO + , N2H+ and 13 CO spectra in clumps GLM1 
(Stage 1), GLM2 (Stage 2), GLM3 (Stage 2), and GLM4 (Stage 2, slightly overlapping with Stage 4 "diffuse red clump"). 
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Fig. 23.— G028.53-00.25: Left: GLIMPSE three-color image, red is 8 /im, green is 4.5 /im and blue is 3.6 /im with MAMBO 1.2 mm 
contours overlaid. The contours are on a log scale from 30 to 320 mjy beam -1 . Right: MIPSGAL 24 /im image with BGPS 1.1 mm 
contours. The contours are on a log scale from 0.12 to 1.3 Jy beam -1 . 
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Fig. 24. — G028. 53-00. 25: Top left: BGPS 1.1 mm continuum dust emission overlaid with linear HCO + contours from 0.43 to 7 K 
km s . Top right: BGPS 1.1 mm continuum dust emission overlaid with linear N2H + contours from 0.9 to 3.6 K km s _1 Bottom left: 
GLIMPSE 8 fim overplotted with red BGPS beam-sized apertures that were used to determine clump masses. Bottom right: HCO + , 
N2H+ and 13 CO spectra in clump GLM1 (Stage 2). The HCO + and N2H+ spectra have been multiplied by a factor of ten, in order to see 
them more clearly relative to the 13 CO. 
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Fig. 25.— G034. 43+00. 24: Left: GLIMPSE three-color image, red is 8 /im, green is 4.5 /im and blue is 3.6 /an with MAMBO 1.2 mm 
contours overlaid. The contours arc on a log scale from 30 to 2400 mjy beam - 1 . Right: MIPSGAL 24 fim image with BGPS 1.1 mm 
contours. The contours are on a log scale from 0.12 to 4.3 Jy beam -1 . 

was observed on the Submillimeter Array by iRathborne et al.l ([20081 ) who found that GLM1 contains at least three 
compact (less than ~ 0.06 pc) cores. This quiescent cloud could be the host of low mass star formation, or it could be 
a transient density enhancement or a quiescent proto-cluster. 

G03443+00.24 

G034. 43+00. 24 (see Figures [231 and |2"6"|) is a stunni ng; example of an active IRDC. This object is just northeast of the 
enormously bright H II region at G34.3+0.2, RH85 ([Reid fc Ho 1985). This complex contains the bright H II region, 
RH85, and the shells of supernova remnants. G034. 43+00. 24, along with a few other IRDCs not discussed in the paper, 
are extremely filamentary structures, streaming radially out of this bright H II region complex, increasing in width 
with distan ce from the complex . This is an excellent source in which to study the triggered formation mechanism of 
IRDCs fe.g lRedman et al.ll2003l) . 

South of G034.43+00.24: GLM2 is long train of "green fuzzies," along the IRDC. IRathborne et al.l ()2008D observed 
G034. 43+00. 24: GLM1 with the Submillimeter Array and found that GLM1 remains unresolved and contains a core 
smaller than the ~ 0.03 pc beam. GLM2 itself is an excellent example of a "diffuse red clump," offset from the 
millimeter peak. The bright 8 fim emission nullifies the use of an extinction mass in this clump; the extinction mass 
is 60 M© while the BGPS 1.1 mm mass is 780 M©. This clump is one of the most "massive" sources in our sample 
according to the BGPS 1.1 mm mass, which is almost certainly biased by a higher dust temperature than the assumed 
15 K. This spectacular clump is an example where all of our conventional mass tracers fail. The HCO + and N 2 H + arc 
extremely bright in all the clumps, especially GLM1 and GLM2 (the brightest by far in our sample). There is also 
strong self- absorption in 13 CO toward these clumps, further indication of the high density and temperatures of these 
active clumps. 
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Fig. 26. — G034. 43+00. 24: Top left: BGPS 1.1 mm continuum dust emission overlaid with linear HCO + contours from 0.55 to 38 
K km s^ 1 . Top right: BGPS 1.1 mm continuum dust emission overlaid with linear N2H+ contours from 0.48 to 27 K km s —1 Middle: 
GLIMPSE 8 fim overplotted with the BGPS beam-sized apertures that were used to determine clump masses. The blue X's are VLA 3.6 
cm point sources, and extended 3.6 cm sources are depicted as blue contours in the adjacent box. Right: HCO+, N2H+ and 13 CO spectra 
in clumps GLM1 (Stage 3), GLM2 (Stage 4), and GLM3 (Stage 2). 



